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ABSTRACT 

We model the production of OH"*", H2O"'", and HaO^" in interstellar clouds, using 
a steady state photodissociation region code that treats the freeze-out of gas species, 
grain surface chemistry, and desorption of ices from grains. The code includes PAHs, 
which have important effects on the chemistry. All three ions generally have two peaks 
in abundance as a function of depth into the cloud, one at Ay ^ 1 and one at Ay ~ 
3 — 8, the exact values depending on the ratio of incident ultraviolet flux to gas density. 
For relatively low values of the incident far ultraviolet flux on the cloud {x ^ 1000; 
X = 1= local interstellar value), the columns of OH^ and H20''~ scale roughly as the 
cosmic ray primary ionization rate C,crp divided by the hydrogen nucleus density n. 
The HsO^ column is dominated by the second peak, and we show that if PAHs are 
present, A(H30'^)~ 4 x 10^'^ cm^^ independent of Ccrp or n. If there are no PAHs or 
very small grains at the second peak, A^(H30"'") can attain such columns only if low 
ionization potential metals are heavily depleted. We also model diffuse and translucent 
clouds in the interstellar medium, and show how observations of A(OH"'")/A^(H) and 
A'(0H"'")/A'"(H20"'") can be used to estimate Qcrp/n, xl'^ and Ay in them. We compare 
our models to Herschel observations of these two ions, and estimate Ccrp~ 4 — 6 x 
10"i6(n/100 cm-3) g-^ and x/n = 0.03 cm'^ for diffuse foreground clouds towards 
W49N. 

Subject headings: astrochemistry, ISM: cosmic rays, ISM: molecules, ISM: clouds, galax- 
ies: ISM, submillimeter 



1. INTRODUCTION 



The OH"'", H2O"'", and HsO"'" ions form the backbone of interstellar chemistry and are important 
probes of the cosmic ray ionization rates in diffuse clouds, on the surfaces (Ay < 1) of molecular 
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clouds, and (with less reliability as we will discuss in this paper) in the interiors {Ay ~ 4 — 8) of 
molecular clouds. They are the backbone of chemistry because once H2 forms, cosmic ray ionization 
of H or H2 leads to the formation of these ions, and HsO^ recombination with electrons leads to 
OH and H2O. Reaction of OH with C or C+ leads to CoQ Once these basic molecules are formed, 
many of the other polyatomic and rare species follow. 

The pathway of cosmic ray ionization of hydrogen to these molecular ions follows two routes 
(see Figure 1). In gas with significant H atoms the ionization of H leads to H"^ that then proceeds 
via a series of reactions (see Figure 1 top) to 0H+, H20^ and HaO"^ ions. We note that the charge 
exchange of O with H^ is slightly endothermic, so the reaction rate is proportional to exp(-230 
K/T); this means that this reaction slows down at cooler temperatures, and a greater fraction of the 
cosmic ray ionizations of H are followed by recombination of H"*" with neutral or negatively charged 
polycyclic aromatic hydrocarbons (PAHs or PAH~) or electrons rather than proceeding to form O"*" 
and then 0H+ . It is this atomic route to OH"'' which is primarily important in diffuse clouds and in 
the Ay ^ 2 surfaces of molecular clouds. A second route dominates deeper in the opaque interiors 
of molecular clouds. Here the ionization of H2 leads to H^, which then proceeds via a series of 
reactions (see Figure 1 bottom) to OH^, H2O''' and HsO"*". A key competitor here to the formation 
of 0H+ is the dissociative recombination of Hj^ with electrons, which has a high rate coefficient 
compared to the reaction of H^ with O. In addition, the reaction of Hj^ with CO dominates that 
with O when the CO abundance exceeds that of O. Therefore, low electron abundances and high 
O abundances are needed to ensure that a large fraction of cosmic ray ionizations of H2 eventually 
produces OH"*". 

As we shall show in this paper, these two routes lead generally to two peaks in the abundances of 
0H+, H2O+, and H3O+ as a function of depth or Ay into a cloud. The first peak (at Ay ~ 0.01 — 1 
depending on the ratio of the incident FUV flux to the gas density) occurs in atomic gas where 
the cosmic ray ionization of H begins the chemical chain|§ The second peak (at Ay ~ 3 — 8, 
again dependent on the FUV flux/gas density ratio) occurs in molecular gas where the cosmic 
ray ionization of H2 begins the chemical chain. Deeper in the cloud the gas phase oxygen freezes 
out as water ice on grain surfaces (e.g., Hollenbach et al 2009, hereafter H09), and the gas phase 
abundances of the three ions drop. Whereas significant columns of 0H+ and H20'^ are produced 
in the first peak, most of the HsO^ column arises in the second peak. 

In warm (T > 300 K) neutral gas with significant abundances of H2 there are other dominant 
pathways to form 0H+, H20^, and HsO^ ions than the two paths initiated by cosmic rays and 
shown in Figure 1. One of the goals of this paper is to show under what conditions (e.g., FUV flux, 
gas density, cosmic ray ionization rate) cosmic rays initiate the formation of the OH"^, H2O"*", and 
HsO"*" ions, and under what conditions other chemical routes dominate. 



CO also has another route initiated by the production of the CH ion, which we do not discuss in detail in this 
paper although it is included in our models. 

^FUV is defined as photons in the range 6eV< hv < 13.6 eV. 
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If cosmic rays dominate, then the observed columns of OH"*", H20''~, and HsO^ serve as a 
probe of the cosmic ray ionization rate. The major goal of this paper is to show how these columns 
depend on the cosmic ray ionization rate, the gas density, the FUV flux, the total column or optical 
extinction Ayt through the cloud, and the molecular hydrogen abundance. Comparison with our 
models allow the cosmic ray ionization rate to be estimated from these ion columns and the atomic 
H column. 

Our method of estimating cosmic ray ionization rates from the abundance of eventual products 
of this ionization dates back to the 1970's, when Black & Dalgarno (1973) pointed out the sensitivity 
of the OH and HD abundances to cosmic ray ionization rates. A later comprehensive paper by van 
Dishoeck &: Black (1986) showed how the abundance of OH in diffuse or translucent clouds could 
be used to estimate the cosmic ray rate. One advantage of using OH"*" or H2O"'" to probe cosmic 
ray ionization rates over using OH is that OH is produced after a series of reactions following the 
formation of 0H+ (see Figure 1). This makes the inference of cosmic ray ionization rate dependent 
on knowledge of the rate coefficients of these additional reactions. 

Recent observations strongly motivate the theoretical study of the OH"*", H20^, and HaO"*" 
molecular ions. Prior to 2010, the only observations of these three ions within the interstellar 
mediunjfl consisted of a few detections of HsO"*" in rich interstellar sources of submillimeter line 
emission and absorption (Wootten et al 1991; Phillips, van Dishoeck & Keene 1992, hereafter 
PvDK92; Goicoechea &: Cernicharo 2001; van der Tak et al 2006). Thanks largely to absorption-line 
spectroscopy with the HIFI instrument on the Herschel Space Observatory, together with additional 
ground-based detections of OH"*" obtained with the APEX telescope (submillimeter rotational tran- 
sition) and the ESO Paranal observatory (near-UV electronic transition), the observational picture 
has improved radically over the past two years. All three molecular ions have now been detected in 
both the diffuse and dense Galactic interstellar medium (Gerin et al. 2010; Ossenkopf et al. 2010; 
Wyrowski et al. 2010; Neufeld et al. 2010; Schilke et al. 2010; Gupta et al. 2010; Bruderer et al. 
2010; Benz et al. 2010; Krelowski, Beletsky, &; Galazutdinov 2010) and in external galaxies (van 
der Tak et al 2008; Weiss et al. 2010; van der Werf et al. 2010; Gonzalez- Alfonso et al. 2010; Aalto 
et al 2011). The most reliable column density determinations are obtained for foreground molec- 
ular clouds lying along the sight-lines to bright submillimeter continuum sources in the Galactic 
disk. Here, the absorbing material typically covers a wide range of line-of-sight velocities, arising 
in widely distributed material along the sight-line. The total inferred column densities of the ions 
lie in the range few x 10^^ — few x lO^^cm"-^, the largest values being attained for OH"*" and the 
smallest for H3O+ (Gerin et al. 2010; Neufeld et al. 2010). 

In their study of OH^ and H2O''' absorption along the sight-line to the luminous star-forming 
region W49N, Neufeld et al. (2010, hereafter NIO) measured an average 0H"'"/H20^ abundance 
ratio of ~ 10, with variations over the range ~ 3 — 15. These observed ratios are considerably larger 



^Note, however, that the H2O''' ion has long been observed in cometary comae (e.g. Wehinger et al. 1974), where 
it is produced by the photoionization of water by solar ultraviolet radiation. 
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than the value ~ 1 expected in fuhy molecular gasO By means of a simple analytical treatment 
of the chemistry of OH"'" and H2O"'", confirmed by more detailed pure gas-phase models performed 
using the Meudon PDR model (Le Petit et al 2006; Goicoechea & Le Bourlot 2007), NIO concluded 
that the molecular fraction in the absorbing material lies in the range 2 - 8 %. This conclusion was 
supported by the observed distribution in velocity space of the OH"'" and H2O"*' absorption, which 
proved similar to that of the atomic gas probed by 21 cm absorption studies, and quite dissimilar 
from that of the molecular gas traced by HF or CH absorption. The analytic treatment introduced 
by NIO also allowed the cosmic ray ionization rate to be inferred from the abundances of OH"'" 
and H20^ relative to atomic hydrogen. The resulting estimate of the cosmic ray ionization rate in 
the range Ccrp= 0.6 - 2.4 xlO~^^ s~^ (primary ionization rate per H atom) was broadly consistent 
with earlier values inferred quite independently from observations of the H3" molecular ion toward 
different sight lines (Indriolo et al. 2007). One of the goals of the present study is to refine the NIO 
analytic treatment of diffuse cloud chemistry through detailed modeling, as an aid to interpreting 
the growing body of observational data now available. 

This paper is organized as follows. In §2 we describe the chemical/thermal model of both 
an opaque molecular cloud illuminated by FUV radiation, or a diffuse cloud illuminated by the 
interstellar radiation field. In §3 we show the model columns and column ratios of OH"'", H2O"'", 
and H3O"'" as functions of the cloud gas density n, the incident FUV flux x (in units of the local 
average interstellar field, see below), and the primary cosmic ray ionization rate C,crp per H atom 
for the case of opaque molecular clouds. We also show the same results for diffuse clouds, but 
with the total column or Ayt through the cloud as an additional parameter. In §4 we compare our 
model results with previous H3O"'" observations and recent 0H+, H2O"'", and H3O"'" observations by 
Herschel, and show how observations compared to models can constrain the cosmic ray ionization 
rates. We summarize our results in §5. Appendix A presents tables of key reaction rate coefficients 
and adopted abundances and grain parameters. Appendix B includes analytic expressions that 
explain the variation of OH"*" with Ay over the first peak, and the variation of H3O"'" with Ay 
over the second peak. Appendix C assesses the sensitivity of our results to certain chemical rate 
coefficients, the gas phase abundance of elemental oxygen and low ionization potential metals, and 
the freeze-out of species. 



^In fully molecular gas the formation rates of these ions per unit volume are equal because every OH^ formation 
is followed by H20^ formation and the destruction rate per ion is nearly the same (due to reaction with H2), leading 
to similar abundances for both ions (NIO). 
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2. THE CHEMICAL AND THERMAL MODEL OF A CLOUD 

2.1. Summary of Prior PDR Model and Modifications 

The numerical code we have developed to model the chemical and thermal structure of an 
opaque cloud externally illuminated by FUV flux is based on our previous PDR model described in 
H09. This ID code models a constant density n (the hydrogen nucleus density) slab of gas, illumi- 
nated from one side by an FUV flux of 2.7 x 10~^x srg cm~^ s~^ incident perpendicular to the slab. 
The unitless parameter x is defined above in such a way that Y ~ 1 corresponds to the average 
local interstellar radiation field in the FUV band (Draine 1978) O The code calculates the steady 
state chemical abundances and the gas temperature from thermal balance as a function of depth 
into the cloud. It incorporates 63 chemical species, ~ 300 chemical reactions, and a large number 
of heating mechanisms and cooling processes. The chemical reactions include FUV photoionization 
and photodissociation; cosmic ray ionization; neutral-neutral, ion-neutral, and electronic recombi- 
nation reactions. H2 self-shielding is included as described in HOD, and CO self-shielding and the 
partial shielding of CO by H2 is included as described in Visser et al (2009). The code includes 
the photodissociation of molecules by "secondary" FUV photons produced (ultimately) by cosmic 
rays (Prasad &: Tarafdar 1983). We also include reactions with charged dust grains and PAHs; 
and the formation of H2, OH, H2O, CH, CH2, CH3, and CH4 on grain surfaces. The code treats 
the freezing of all condensable species to grain surfaces and three desorption processes: thermal 
desorption, photodesorption, and cosmic ray desorption. The only significant difference in the 
desorption code used here versus the H09 code is the inclusion of the new (higher by factor ~ 3) 
rate of photodesorption of CO (Oberg et al 2009). The code does not include photodesorption 
by the secondary FUV photons; this process is negligible at the peaks in the 0H+, H2O"'', and 
HsO"*" abundances. The code has been used to model regions which lie at hydrogen nucleus column 
densities iV < 4 x 10^^ cm~^ (or Ay ^ 20) from the surface of a cloud. Therefore, it applies not only 
to the photodissociated surface region, where gas phase hydrogen and oxygen are nearly entirely 
atomic and where gas phase carbon is mostly C"^, but also to regions deeper into the molecular 
cloud where hydrogen is in H2 molecules and carbon is in CO molecules. Even in these molecular 
regions, the attenuated FUV field can play a significant role in photodissociating H2O and O2, in 
photodesorbing species adsorbed on grain surfaces, and in heating the gas. However, the code is 
now sufficiently general that it finds the steady state solutions for abundances and temperature in 
any region of a molecular cloud, even where FUV is insignificant. 

We emphasize that we present a steady-state model of chemical abundances as a function of 
depth into the cloud. The chemical timescales can be quite long, which might suggest that time 
dependent models are more appropriate. For example, the timescale to convert atomic H gas to fully 



^ Note that x — 1 corresponds to Go = 1.7 in the Tielens & HoUenbach (1985) units based on the Habing (1968) 
local interstellar radiation field. The shape of the FUV spectrum for x > 1 is implicitly assumed to mimic that of 
the Draine field, which is approximately that of a Te// ~ 30, 000 K star. 
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molecular H2 gas is years, where n is the hydrogen nucleus density in units of cm . 

However, as we show below, the OH"*" and H2O"'" ions peak in abundance when a:(H2) ~ 0.03 — 0.1, 
which occur typically at Ay ~ 0.1. [Note, abundances in this paper are defined relative to hydrogen 
nuclei, so that x(H2) = n(H2)/n and n ~ n(H) + 2n(H2)]. Therefore, the timescales for even diffuse 
clouds of density n ~ 100 cm~^ to reach these abundances are less than ~ 10^ years, which is shorter 
than the typical lifetime of a diffuse cloud (Wolfire et al 2003). Steady state solutions therefore 
generally apply, at least for computing the columns of these ions. The steady state models may 
somewhat underestimate the ion columns for low density diffuse clouds with Ay > 0.1, since the 
steady state solutions can lead to lower abundances of these two ions than time dependent solutions 
in the high Ay regions. This arises because the steady state abundances of H2 are higher than 
time dependent models which start with fully atomic gas. Higher H2 abundances lead to more 
destruction of OH"'" and H2O"'" in the high Ay regions where H2 and not electrons dominate the 
destruction of these ions. Liszt (2007) provides a detailed analysis of the time dependent formation 
of H2 and OH"*" as a function of Ay in diffuse clouds with initial atomic conditions. 

We also emphasize that our model does not include turbulent dissipation and heating of small 
pockets of gas along the line of sight (e.g., Godard, Falgarone, & Pineau des Forets 2009 and 
references therein). We do, however, run PDR models with small fractions of the line of sight 
having either enhanced temperature or enhanced rates of ion-neutral drift to test the possible 
effects of turbulence, and find the effects on 0H+, H20^, HsO^ column densities are likely small. 

For simplicity we assume constant H nucleus density n in our models. At low Ay and low 
a;(H2) < 0.1, where much of the 0H+ and H2O'*' columns often arise, the temperature is quite 
constant so that constant density implies constant thermal pressure. If thermal pressure (and not 
turbulence) dominates deeper into a high Ay cloud, then the density will rise as one moves from 
its warmer PDR surface to its CO-cooled molecular interior. In addition, self gravity can raise the 
density of the interior regions at higher Ay . Still another effect is the transition of atomic hydrogen 
to molecular hydrogen which can raise n by a factor of 2 if thermal pressure is conserved. We ignore 
the possible rise in density, which mainly affects the second peaks of the ions deep {Ay > 4) in 
the cloud. If such a density enhancement occurs, it tends to depress the abundances of 0H+ and 
H2O"'" in the second peaks, but the abundance of HsO^ at the second peak is not sensitive to 
hydrogen density if PAHs are present (see §3) If] The second peak is dependent on the abundance 
of PAHs, very small grains and low ionization potential metal ions there, which can control the 
electron density, n^. The abundance of PAHs and very small grains at high Ay is uncertain due to 
their possible coagulation on larger grain surfaces. 

One significant difference between the chemical code used in this paper and that used in H09 
is the inclusion of PAHs and very small grains (VSGs, radius ^50 A), which affect the ionization 
balance by enhancing the recombination of positive atomic ions. In the rest of this paper we shall 



^If a reader is interested in the second peak and knows the density n of a given source there, then our models with 
that n will give a good prediction of the behavior of the second peaks. 
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often use the term "PAH" to denote both PAHs and VSGs. PAHs also affect the second peaks 
of the ions by destroying electrons in the reaction e + PAH — t- PAH~. (The effects of PAHs on 
ionization balance has been treated extensively before in the literature: e.g., Lepp & Dalgarno 
1988, Bakes k Tielens 1998, Weingartner & Draine 2001a, Flower & Pineau des Forets 2003, Liszt 
2003, Wolfire et al 2003, 2008). H09 focused on the HgO and O2 peaks, which occur at relatively 
high Ay ~ 3 — 7. Although the PAH abundances are not known deep in molecular clouds, H09 
assumed that their abundances are low, due to coagulation of PAHs on the surfaces of larger dust 
grains. Here, however, we are mostly interested in the chemistry at low Ay < 1, and in particular 
in the chemistry of diffuse clouds. These are the regions where PAHs are observed to be present, 
and we use PAH parameters derived from the literature. For simplicity, we adopt a single size PAH 
for the PAH distribution, and assume that the standard PAH has 100 C atoms, and a number 
abundance of xpah = 2 x 10~^ with respect to H nuclei (see Wolfire et al 2008 for a discussion of 
the amount of carbon in PAHs; here we adopt 100 C atoms and not the 35 C atoms that Wolfire 
et al adopted because of the result of Draine & Li 2007, which suggested that the distribution in 
mass peaks at 100 C atoms). One of the key impacts of PAHs is in the recombination of H+. As 
shown in Figure 1 (top) , cosmic ray ionization of H can lead to 0H+ , but a competing route is the 
neutralization of H"*" by an electron, PAH, or PAH~. PAHs therefore lower the production of OH''" 
and the columns of all the ions in the first peak. Conversely, in the second peaks, the reduction 
in electron abundance caused by PAHs cause fewer Hj]" ions to recombine with electrons, and lead 
to greater production rates of the ions as well as smaller destruction rates for HaO^. Therefore, 
the columns of the ions increase with the presence of PAHs in the second peaks. We also discuss 
results with no PAHs or very small grains at high Ay- 

In Appendix A we present Table 1 which lists the rate coefficients of key reactions in the 
pathways to the 0H+, H2O"'", and H3O"'' ions, as well as reactions that are either new or changed 
since H09. Of particular note here are the photoionization of OH and H2O, the photodissociation of 
0H+ and H2O+, the treatment of PAHs-especially the photodetachment reaction rate for PAH~, 
the fine structure level population dependence in the charge exchange reaction of O with H"*", and 
some minor modifications in reaction rates key to determining the abundances of the OH^", H20^, 
and HaO"*" ions, such as their dissociative recombination rates with electrons and their reactions 
rates with H2. 

We have also included in Table 1 rates that are important in producing H^ by chemical means 
rather than by cosmic ray ionization of H in PDRs. There are two main chemical routes. The first 
is initiated by the FUV photoionization of C to C+. The C+ reacts with OH to form C0+. The 
C0+ charge exchanges with H to form H+ . The second is also initiated by FUV photoionization of 
C to C+. The C+ reacts with H2 to form CH+. The CH+ is photodissociated by FUV photons to 
produce H+. Both these reaction chains are very much enhanced by high ( > 300 K) temperatures 
in gas with appreciable H2. The first chain is enhanced because high gas temperatures lead to 
significant amounts of OH being formed by the reaction H2 + O — )■ OH + H. This reaction has 
an activation barrier of AE/k = 3160 K and so is insignificant for low gas temperatures. The 
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second chain similarly is enhanced because the reaction C"*" + H2 ^ CH"*" + H has an activation 
barrier of /^E/k = 4640 K (see Table l)|l| In regions where either of these two chains dominate the 
production of H"*" over cosmic ray ionization, the ions OH"'", H2O"'", and HaO"*" only provide upper 
limits to the cosmic ray ionization rates. 

Other alternate routes to the production of OH"^ and 1120"*" ions include three routes which 
produce OH or H2O without the ions: the production of water on grain surfaces followed by 
photodesorption of the water to produce gas phase OH and H2O, the radiative association of O 
with H to form OH, and the reaction of FUV-pumped H2 in excited vibrational states with O to 
form OH. The gas phase OH and H2O can be photoionized by FUV photons to produce OH"*" and 
H20"^ (see Table 1). Unlike the routes described in the preceeding paragraphs, these routes are 
never dominant in producing the ion column densities. 

An analogous route to the formation of H2O on grain surfaces followed by photodesorption, but 
one not treated in this work, is the time-dependent evaporation of water ice which occurs around 
newly formed stars. The rapid rise in embedded luminosity heats the dust grains above about 100 
K, and the icy mantles on grains are then thermally evaporated. This sudden injection of high 
abundances of water vapor into the gas is followed by reaction of the gas phase water with HCO"*" 
and H|^ to form H3O"*". Eventually, the system relaxes to the steady state chemistry described in 
this paper, but for a short time, there might be a large enhancement in HsO"*" (Millar, Herbst, & 
Charnley 1991, PvDK92). If this release of H2O from the icy grain surface to the gas occurs in 
regions with elevated FUV fields, the photoionization of H2O could also result in enhanced H2O"'" 
abundances (Gupta et al 2010) 

We also present in Appendix A Table 2, which lists the gas phase elemental abundances, the 
PAH properties, and the grain surface area per H nucleus adopted in our code. 

2.2. Interstellar Cloud Models 

As noted above, our opaque molecular cloud models invoke a constant density n slab, illumi- 
nated from one side by a ID normal FUV flux x- We solve for the gas temperature and the gas 
phase and ice abundances of the various species as a function of depth in the slab. Note that depth 
is synonymous with hydrogen nucleus column N [=A^(H)-|-2A^(H2)] into the cloud or Ay into the 
cloud. We take N = 2 x 10^^ Ay cm~^. We follow the chemistry to Ay ~ 20, beyond which there 
is little contribution to the columns of OH"*", H2O"'", and HsO"*" ions. 

We use the primary cosmic ray ionization rate per H atom as an input parameter, since our code 
calculates the secondary ionizations caused by cosmic rays and these depend on the H2 and electron 
abundances. In order to probe the range of cosmic ray ionization rates suggested in the literature. 



^As discussed in Tielens & Hollenbach (1985), our code also includes the reaction of vibrationally excited H2 with 
to form CH+. This reaction has no activation barrier and can be moderately important at low Ay ^0.6. 
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we include cases with primary ionization rates of Ccrp= 2 x 10" s~ and 2 x 10~ s~ per H atom, 
which correspond to total rates (including secondary ionizations) of about d^crt ~ 3 — 5 x 10~^^ s~^ 
and 3 — 5 X 10^^^ s~^ respectively. Note that the primary rate per H2 molecule is 2^crp- We also 
adopt (crp as the primary rate for He atoms; however, He has insignificant secondary ionizations. 
Although there is evidence that the cosmic ray ionization rate may decrease with depth (e.g., 
Rimmer et al 2011, Indriolo Sz McCall 2012, and discussion later in this paper), for simplicity we 
assume the primary cosmic ray ionization rate does not vary with depth into a cloud in a given 
model. The main effect of this assumption is that we may have overestimated the abundances 
and columns of the 0H+ and H2O+ ions in the second (deeper) peak relative to the first peak. 
However, since we vary Qrp in our parameter study, the reader can use lower Qrp for the columns 
in the second peak if so desired. The abundance of H3O"'" in the second peak is not sensitive to 
Ccrp if PAHs are present, as we will show below. 

The main parameters that we explore for our one-sided, opaque molecular cloud models are the 
gas density n, the incident FUV flux %, and the primary cosmic ray ionization rate (^crp- We study 
the parameter space 10 cm~^< n < 10^ cm~^, 1< x < 10^, and 2 x 10"^''' s~^ <Ccr39< 2 x 10""^^ 
s~^. We explore the sensitivity of the columns of 0H+, H2O+, and H3O+ ions to assumptions 
about the PAH chemistry, the elemental gas phase abundances of low ionization potential metals, 
and the rate coefficient for the formation of H2 on grain surfaces. 

The diffuse cloud models treat a constant density slab of total thickness Ayt illuminated 
on both sides by a ID normal FUV flux x/2- We explore the parameter space 10^^^'^ s^^ < 
Ccrp/n2 < 10-^^-^ s-^, 0.01 < Avt < 3, 30 cm-^< n < 300 cm-^, K x < 10, and 1 < x/n2 < 10 
where n2 = n/100 cm~^. The models of Wolfire et al (2003) provide a good estimate of x/n2 
in the Galaxy by estimating x from the star formation rate and the dust opacity as a function 
of galactocentric radius R, and setting the thermal pressure to provide two phases, a cold diffuse 
cloud phase and a warm intercloud medium that fills most of the volume. Assuming R = 8.5 kpc 
as the solar location, Wolfire et al find x = 1-0, n = 33 cm~^, and x/^2 = 3.0 for diffuse clouds 
at R = 8.5 kpc; x = 2.35, n = 49 cm~^, and x/"'2 = 4.8 at i? = 5 kpc; x = 3.0, ra = 54 cm~^, 
and = 5.5 at i? = 4 kpc; and x = 3.8, n = 60 cm~^, x/1^2 = 6.3 at i? = 3 kpc. Therefore, 

our prime region for study is x/"2 = 3 — 6. We are especially interested in how the ratio of the 
columns iV(0H+)/iV(H20+) and iV(OH+)/iV(H) vary as functions of Ayt, x/n and C,crp/n. 

We note that the columns our models predict are columns perpendicular to the 1-sided (molec- 
ular cloud) slab or 2-sided (diffuse cloud) slab. If the slabs arc viewed at an angle 9 with respect 
to the normal, the observed columns will increase by (cos^)^^. Another effect that will obviously 
raise the columns is if there are more than one diffuse cloud (in a given velocity range) along the 
line of sight (los), or if the molecular cloud is clumpy and FUV scattering then introduces several 
"surfaces" along the los. However, as we will show, Ccrp/n in diffuse clouds can be estimated from 
our models from the ratios iV(0H"'")/iV(H20"'") and A^(OH"'")/iV(H), and the ratios arc indepen- 
dent of the geometric effects. Nevertheless, there is some degeneracy in the solution for Ccrp/n, 
depending on the combination of the Ayt of a single cloud, x/f^i and the enhancement in columns 
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created by the geometrical effects. 

3. MODEL RESULTS 
3.1. PDR surfaces of opaque molecular clouds 

3.1.1. The chemical and thermal structure of individual clouds 

In order to understand the columns of 0H+, and HsO^ ions produced as a function 

of n, X-, and Ccrp, we first study the detailed chemical and thermal structure of a few specific 
(standard) cases. Figure 2 shows the chemical abundances as a function of depth Ay into the cloud 
for the case n = 10^ cm~^, x = 1; and C,crp= 2 x 10"^'' s~^. This case is chosen not only because it 
may be appropriate for the ambient interstellar radiation field (ISRF) incident on a relatively low 
density GMC, but also because the surface {Ay ^3) structure (i.e., T and chemical abundances as 
a function of depth or Ay) is illustrative of the depth structure of a diffuse or translucent cloud. 
In addition, the cosmic ray rate may be appropriate to the interior of molecular clouds. 

The main chemical result is that the OH^, H2O"'', and HaO^ ions all have a peak at 0.03 < 
Ay < 0.3, and then have a second peak at Ay ~ 6 for this combination of n and x- We first discuss 
the surface peaks at low Ay. 0H+ peaks at Ay ~ 0.03, where the molecular hydrogen abundance 
is x(H2) ~ 0.01. The H2O+ peaks slightly deeper, at Ay ~ 0.1, where x(H2) ~ 0.1. Finally, H3O+ 
peaks at Ay ~ 0.3, where x(H2) ~ 0.25. Appendix B provides an approximate analytic solution 
to the chemistry that explains the 0H+ first peak, and its relation to x(H2). From Ay ~ 0.01 to 
either Ay ~ 0.03 (0H+) or Ay ~ 0.1 (H2O+) or Ay ~ 0.3 (H3O+) the three ions rise in abundance 
with increasing Ay because of the rise in the H2 abundance, which drives the cosmic ray-produced 
O"*" to the molecular ions. At larger Ay, but for x(H2) significantly less than its maximum value of 
0.5, the abundance of 0H+ and H20^ plateau at the peak value because here both their formation 
and destruction rates are proportional to x(H2). H2O''' tends to peak at somewhat higher Ay than 
0H+ because, even at the peak of OH^ abundance, not all cosmic ray ionizations lead to H20^ and 
so its abundance continues to rise as the H2 abundance rises with increasing Ay. Finally, as x(H2) 
approaches 0.5, two effects lead to a drop in the 0H+ and H2O+ abundances with increasing Ay. 
One is that the formation rates of these ions saturate as nearly every cosmic ray ionization leads to 
their production, whereas the destruction rates still scale as x(H2), which increases with increasing 
Ay. The other dominant effect is that the H abundance drops so that the 0H+ formation rate 
via the top chain of reactions in Figure 1 drops. The bottom chain at relatively low Ay is not as 
efficient at producing OH^, as the electrons are relatively abundant (x^ ~ 10"'*) in these surface 
regions, and Hj^ recombines with electrons rather than forming OH^|j The HaO"*" is destroyed 



*As a result of the inefficiency of the bottom chain when the electron abundance is relatively high, the first peak 
is always dominated by the top chain of reactions in Figure 1, even as x{Y{2) approaches 0.5. The bottom chain 
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not by H2, but by electrons, whose abundance stays quite constant (supplied mostly by C"*" but 
with possible contribution by H+ at high Ccrp/n) at the cloud surface. The HsO"*" starts to drop 
in abundance once the gas becomes predominantly H2, due to the second effect described above. 

We next discuss the second deeper peak in the OH^, H2O''', and HsO^ ions. As one moves 
deeper into the cloud (typically, Ay ^2), the electron abundance starts to drop. As this happens, 
a greater fraction of the cosmic ray ionizations of H2 leads to the production of the three ions, and 
their formation rates rise. The destruction of 0H+ and H2O+ is by H2, which now has constant 
abundance (the gas is fully H2), so the destruction rates hold constant. Therefore, these two 
ions rise in abundance. Finally, they peak and fall in abundance for Ay > 6 because gas phase 
oxygen freezes out as water ice, and again the oxygen reaction with Hjj" cannot compete with Hg" 
electronic recombination or its reaction with CO. Therefore, the formation rates of all three ions 
drop. H3O+ behaves somewhat more dramatically, because its destruction is mainly by dissociative 
recombination with electrons. Thus, as the electron abundance drops, not only is the formation 
rate of HsO"^ enhanced, but the destruction rate is suppressed. Therefore, H3O"'" rises to much 
higher abundances than 0H+ and H2O"'" in the second peak. It also drops at very high Ay because 
gas phase elemental oxygen from which H3O"'" forms freezes out as water ice. 

In this particular case, there is a column of A/"(0H+)=2 x 10^-*^ cm~^ in the first surface peak 
and 4 x 10^^ cm~^ in the deeper peak; for H2O"'" the columns are 1.4 x 10^^ cm~^ and 2 x 10^^ 
cm^^; and for HaO^ the columns are 1.1 x 10^^ cm~^ and 9.4 x 10^'^ cm^^. The columns of OH"*" 
and H20^ are not detectable. Typically, columns of > 10^^ cm~^ arc needed for detection via 
absorption spectroscopy. In particular, observations of 011+ often imply columns > 10^^ cm~^, 
which suggests that higher cosmic ray rates are required. Therefore, for the rest of our standard 
cases we use Ccrp= 2 x 10~^^ s~^. Most of the H3O+ column is produced in the second deeper 
peak, and in our model the predicted columns are detectable. For an absorption measurement a 
background submillimeter source behind or in a cloud with Ay > 6 from observer to submillimeter 
source is required. 

Figure 3 shows the case n = 10^ cm~^, x = 1) ^-nd Ccrp= 2 x 10""*^^ s~^: in other words, 
the same as Figure 2 but with ten times the cosmic ray ionization rate. The abundances of all 
three ions in the first peak rise in proportion to (crp ■ The 0H+ and H2O''" ion abundances in the 
second peak also scale roughly with (crp- The column of 7V(OH+)=2.2 x 10^^ cm~^ in the first 
peak and 3.2 x 10^^ cm~^ in the deeper peak; for 1120"*" the columns are 1.5 x 10^^ cm~^ and 
6.9 X 10^^ cm~^; and for HsO"*" the columns are 9.2 x 10^^ cm~^ and 7.3 x lO^'^ cm~^. Because 
the electron abundance is low in the second peak, a significant fraction (.1 — 0.3) of cosmic ray 
ionizations of H2 lead to 0H+ and H2O+, and their destruction is by H2, which does not change 
in abundance with varying (crp- Thus, the scaling of 011+ and 1120+ abundances and columns 



dominates in the second peak, however, because the H2 abundance is high and the electron abundance is low. 
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with CcrpU However, the HsO^ abundance in the second peak does not rise hnearly with Ccrp, but 
stays fairly constant, because although the formation rate scales with C,crp, the destruction rate also 
increases as Ccrp is raised, due to the higher electron abundances produced by the enhanced cosmic 
ray flux. Unlike 0H+ and H2O"'", which are destroyed by H2, H3O+ is destroyed by dissociative 
recombinations with electrons. In fact, to first order, we would expect the abundance of HsO^ to 
scale as Ccrp/ne =C,cTp/ {xen) at the second peak. 

The abundance of electrons Xe deep in the cloud depends on the uncertain PAH abundance 
deep in the cloud. If the PAH abundance remains as high as is indicated in diffuse clouds and cloud 
surfaces (as we assume in our standard models), then we obtain the following result. Electrons are 
formed by cosmic ray ionization of H2. Electrons are mainly destroyed by collisional attachment to 
neutral PAHs, and the PAHs are mostly neutral. Therefore, the abundance of electrons Xe (xC,crp/n. 
As a result, if PAHs are abundant in this deep peak, we predict that x(H30+) oc C,crp/{xen) will be 
independent of both n and Ccrp! Note that Figure 3 compared to Figure 2 shows that the electron 
abundance in the second peak does scale roughly as Ccrp, and that the abundance of H3O+ in 
the second peak does not change significantly as we increase the cosmic ray ionization rate by 10. 
Appendix B presents an analytic solution for ^(HaO"'") near the second peak if PAHs are present. 
Since cosmic ray ionization of H2 is similar to X-ray ionization of H2, this result implies that, if 
PAHs are present, regions of enhanced X-ray ionization will not show enhanced HaO"*" columns. 
We discuss in §3.1.3 the case of no PAHs at high Ay. 

Figure 4 shows the same case as Figure 3, but plots a parameter e, first discussed by NIO. 
Here, the parameter e is defined as the rate per unit volume of formation of OH^ divided by the 
total (not primary) rate per unit volume of cosmic ray ionization of H and H2. In effect, e is an 
efficiency parameter in determining the formation of OH"^ from cosmic rays. If PAH, PAH~ and 
electron abundances are relatively low, and H2 and gas phase O abundances are high, then e is 
near unity. Essentially, one needs O"*" to react with H2 before H"*" reacts with e, PAH, or PAH~. 
Although this qualitative limit is clear from the chemical pathways shown in Figure 1, we derive 
in Appendix B an analytic formula for e which makes this statement more quantitative. Roughly, 
the condition for e to be of order unity is: 



However, if the reverse is true, then H+ can recombine with PAH, PAH^, or electrons and disrupt 
the chain of reactions that lead to OH"*", leading to low (< 1) values of e. Figure 4 also plots the 
temperature and repeats the plots of the abundances of H2 and electrons, since they help determine 
the value of e, as well as the abundance of OH"*" (see Appendix B). Finally, we add the abundance 
of Hg" to Figure 4 since it is also used to estimate cosmic ray rates (e.g., Indriolo et al 2007, 2011). 



'Even at the second peak, both electrons and CO compete with O in reacting with HJ. 
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Because Hg" is formed by the reaction of H2 with and often destroyed by electrons (see Figure 
1), we see the Hj^ abundance rise with Ay as the H2 abundance rises and the electron abundance 
falls. In general, the probes the cosmic ray ionization rates at higher Ay than the first peak in 
0H+. 

Figures 5 and 6 show the case n = 10^ cm~^, x = 100, and Ccrp= 2 x 10^^® s^^. This case 
is nearly identical to the standard case in H09, and is representative of a Giant Molecular Cloud 
(GMC) surface illuminated by an FUV field somewhat higher than the interstellar radiation field 
because of the presence of nearby O and B stars. The total cosmic ray ionization rate (~ 3— 5 x 10~^^ 
s~^ including secondary ionizations) may represent values on the surfaces of GMCs, but may be 
somewhat high for the interior. We see from Figures 4 and 6 that the gas temperatures of the 
molecular interiors of these clouds at Ay > 5 is of order 30 K for n = 10*^ cm~^ and 70 K for 
n = 100 cm^'^, due to cosmic ray heating when Ccrp= 2 x 10~^^ s~^. Typically, temperatures in 
molecular cloud interiors are observed to be ^ 30 K, suggesting that such high cosmic ray ionization 
rates may not be appropriate for molecular cloud interiors. However, as we shall see, such high 
cosmic ray rates are required to explain observations of diffuse clouds, which should have the same 
cosmic ray rates as the surfaces (Ay ^ 2) of GMCs. This suggests that <^crp may be higher on the 
surface of a molecular cloud than deep in its interior. 

If chemistry is driven by FUV photoreactions and particle-particle reactions, such as H/H2 
chemistry, then the chemical abundances mainly depend on the ratio x/'^- Therefore, one expects 
and sees that the H2 abundance of Figure 5 closely matches that of Figures 2 and 3, which have 
the same x/i^ ratio. The H2 abundance in the higher x case is a bit lower because of FUV 
photodissociation of H2 in FUV-pumped excited vibrational states of H2. 

However, the molecular ion abundances in the first peak depend to first order on the ratio 
Ccrp/n (see Appendix B) and are not too sensitive to %. Thus, in this n = 10^ cm~^ case, their 
abundances in the first peak drop by a factor of nearly 10-30 compared to Figure 3 (which has 
the same value of Ccrp) ^ the density rises by 100 from the n = 10^ cm^'^ assumed in Figure 3. 
Slight differences in electron abundances, H2 abundances, and T explain the divergence from the 
expected dependence. 

The second peaks of OH"*' and H2O''' nearly scale expected. However, the second 

peak of the H3O"'" abundance (~ 10~^) is independent of n, as predicted above by the scalings of 
the electron density with n and Ccrp if PAHs are present. 

Figures 7 shows the case n = 10® cm~^, x = 10^, and Ccrp= 2 x 10^^® s^-*^. This case may 
represent strongly illuminated PDRs such as may occur around embedded compact or ultracompact 
HII regions, or possibly embedded protostars illuminating the opaque walls of outfiow cones. This 
high density and high FUV flux case was chosen because most of the column of all three ions is 
produced not by cosmic ray ionization, but by other chemical reactions. In Figure 7 we see an 
enormous enhancement of 0H+ abundance at Ay 1.6. Here, T ~ 1000 K and at the same time 
the abundance of H2 is moderately high, ^ W~^. At these elevated temperatures, as discussed in 
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§2.1, the H2 can react rapidly with O to form OH or with C"^ to form CH+, leading to reaction 
chains that make 0H+, H2O''", and HsO"*". One of the key heating mechanisms providing this high 
T is the FUV pumping and the H2 formation pumping of excited vibrational levels of H2 , followed 
by collisional de-excitation of these levels which leads to gas heating (e.g., Tielens k. Hollenbach 
1985). The essential point is that the 0H+, H2O+, and H3O+ columns are not provided by cosmic 
ray ionization, and therefore cannot diagnose the cosmic ray ionization rate, except to give an 
upper limit. 



3.1.2. Contour plots of integrated columns of ions 

Figure 8 shows the contours of the integrated (from Ay = to Ay = 20) columns of OH"*" 
ions for a primary cosmic ray ionization rate of Qrp='^ x 10^^^ s~^ per H atom, and plotted as 
functions of n and x- Note that the 0H+ columns include both the first and second peaks. We 
emphasize that these are columns perpendicular to the face of the PDR slab. If clouds are observed 
obliquely, proportionately more column will be in the line of sight. Similarly, if we are observing 
a slab illuminated on both sides, the columns will be raised by a factor of 2 if the slab is quite 
optically thick (Ay » 1) so that first and second peaks occur on both sides. The upper right 
hand portion of the figure is blacked out because radiation pressure, photoelectric emission, and 
photodcsorption forces on dust grains, when x/'n.2 300, drives dust rapidly through the PDR 
(Weingartner &: Draine 2001b). Such high ratios of x/n rarely occur in nature, and require a much 
more detailed PDR code. 

Figure 8 shows two main features. First, at low x ^ 1000 or low n 10^ cm~^, the columns of 
0H+ are roughly proportional to and almost independent of %, as discussed above. There is a 
weak dependence on x for x ^ 1000 because higher x leads to higher gas T, and the 0+ abundance 
rises with exp(-230 K/T). From this figure it is clear that since OH"*" columns of at least 10^^ cm^^ 
are needed to make absorption observations feasible, low density (n<300 cm^'^) clouds should be 
targeted if Ccrp= 2 x 10~^^ s~^. It should also be noted, as is obvious in Figures 2, 3 and 5, that 
both the first and second peaks contribute to the 0H+ column so that there is substantial column 
at both Ay < 1 and Ay > 1. The 0H+ in this case is quite cold, T < 200 K. Second, these relations 
completely break down in the upper right portion (n > 10^ cm~^ and x <; 10^) of this contour plot. 
Here, as discussed above, the high x and n (high density brings the H/H2 interface closer to the 
surface, where the gas is warm) lead to very warm (T ~ 1000 K) H2 near the surface (Ay <^ 2) of 
the cloud. This warm H2 drives reactions which produce OH^ without the need for cosmic ray 
ionization. In this way, observable columns of 0H+ can be formed at these elevated values of x 
and n, and the columns occur at moderate Ay ~ 2. The 0H+ in this case is quite warm, T > 300 
K, and its column is independent of (^crp^ but depends on n and x- 

Figure 9 shows the contours of the integrated (from Ay = to Ay = 20) columns of H2O+ 
ions for a primary cosmic ray ionization rate of Ccr-p=2 x 10~^^ s~^ per H atom. This figure can 
be compared with Figure 8, which treats the same cosmic ray case for 0H+. We focus on regions 
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where these ions might be detectable, that is, for columns > 10^^ cm~^, and in the cosmic ray 
dominated zones of n and x- In these low density regions, comparison of Figure 9 with 8 reveals 
that the column ratios arc of order unity, as might be expected. The formation rates of these two 
ions are nearly the same (stemming from the cosmic ray ionization rate) and the destruction rates 
of these two molecules are nearly the same (via H2 in the regions where most of the columns are 
generated). Note that in Figures 2, 3, and 5 the local 0H+ abundance peaks at lower values of Ay 
than the H2O''' abundance, and that Figures 8 and 9 present integrated columns through a high 
Ay slab. Therefore, we can obtain higher column ratios of these two molecular ions if we truncate 
our slabs to small diffuse clouds with Ay < 1 (see below, §3.2.2). 

We do not provide a contour plot similar to Figures 8 and 9 but with different (^crp because the 
results are simple to describe. For x ~ 1000 the OH^ and H20^ columns scale as Ccrp- For x ~ 1000 
the columns are independent of Ccrp because cosmic rays do not produce the 0H+ and H2O"'". 
Forcing a fit with the expected (to first order) Ccrp/n dependence in the cosmic ray-dominated 
region % < 100, we find roughly (factor of 2) the simple expressions 

and 

.,H.O^,.,„.3,..(12^)(_^^)... ,3, 

We emphasize that these columns are the total 0H+ and H2O+ columns in both the first and 
second peaks. 

We do not present integrated (from Ay = to Ay = 20) columns of H3O+ ions because, as 
discussed above, most of their column arises at high Ay ~ 6 and this second peak column is quite 
constant as a function of n, Xi and Ccrp if PAHs are present at the same abundances as in diffuse 
clouds. With PAHs present in the second peak, we find that the HaO"*" column is slightly dependent 
on the gas density, varying from 10^^ cm~^ at n ~ 100 cm~^ to 10^^-^ cm~^ at n ~ 10^ cm~^. 
At densities of n = 100 cm~^ and with x = 1 and Ccr-p=2 x 10~^^ s~^, we remind the reader (see 
discussion of Figure 3) that Ar(H30+)= 9.2 x 10^^ cm-^ in the first peak; like 0H+ and H2O+ the 
column of H3O"'" in the first peak scales with Ccrp/n so that even lower densities or higher Ccrp are 
needed to make HaO"*" detectable in the first peak (e.g., in a diffuse foreground cloud). 

3.1.3. Sensitivity to PAH abundance and other parameters 

PAHs may not have a high abundance deep in the high Ay regions of a molecular cloud, 
perhaps due to coagulation of PAHs on larger grains (see, for example, discussion in H09). We 
have run a number of cases with the PAH abundance set to zero to see the effect on the second 
peaks. In these runs, we also reduce the number of small grains by setting the minimum grain size 
to 100 A, but we still refer below to this case as the "no PAH" case for high Ay cloud interiors . 
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The results depend on the gas phase abundance of metals, like S, Si, Mg, and Fe; see Table 2 

of Appendix A for our assumed gas phase abundances of these species in the interstellar radiation 
field. Our code follows the freezeout of these species as a function of depth into the cloud, as pho- 
todesorption no longer is able to keep the metal atoms off the grain surfaces. The photodesorption 
yields of these species are not known. We find that for yields greater than 10~^, there are still 
sufficient gas phase metals present at the second peak to supply electrons and suppress the H30''" 
abundance by factors of > 10 relative to the case with PAHs. However, if the photodesorption yields 
are very low so that the gas phase abundances of these metals at the second peak are < 3 x 10"^, 
then we obtain H3O+ columns in the second peak comparable to the case with PAHs. For example, 
assuming n = 10"^ and x = 100; we obtain iV(H30+)= 4 x 10^^ cm^^ if C,crp/n2 = 2 x 10~^® s~^ per 
H atom, and N{Yi.^O'^)= 8 x 10^^ cm~^ if C,crp/n2 = 2 x 10~^^ s~^ per H atom. The main difference 
between the PAH and no PAH case lies in the electron loss mechanisms: with PAHs the electrons 
are lost by encountering PAHs; with no PAHs the electrons are lost by recombining with ions. The 
atomic ions have slower rates of recombination than the molecular ions, so the presence of metal 
atomic ions leads to higher electron abundances and lower H30^ columns. The 0H+ and H2O"'" 
abundances and columns also decrease in the second peak, but by smaller factors (^3), because 
their destruction is not by electrons but by H2, whose abundance does not change. 

In Appendix C we discuss the lack of sensitivity of the results on certain key reaction rate 
coefficients, the assumed abundance of gas phase oxygen, and on the rate coefficient for H2 formation 
on grain surfaces. We discuss the effect of raising the gas phase abundances of metals. We also 
treat the increase in the HaO"*" columns in the second peak if elemental O does not freeze (as water 
ice) at high Ay. 

3.2. Diffuse and Translucent Clouds 

Diffuse clouds have relatively small columns. Ay ^ 1, whereas translucent clouds have columns 
intermediate between diffuse clouds and GMCs, 1 < Ay ^ 5. Although GMCs have sufficient column 
to incorporate both peaks of the OH"*", H2O"'", and H3O+ abundances, diffuse and translucent clouds 
typically (with the possible exception of the highest Ay translucent clouds) only contain the surface 
peak. In fact, a diffuse cloud may have such small Ay that the cloud may truncate the full peak 
of one or all of the ions. Therefore, the parameter Ayt, the total Ay through the cloud becomes 
an important new parameter in determining, for example, ratios of the columns of the ions. 

Both diffuse and translucent clouds have insufficient column to be gravitationally bound, so 
generally they are in thermal pressure equilibrium with the interstellar medium. The local typical 
value of the thermal pressure is nT ~ 3 — 4 x 10^ cni^"^ K, and this pressure may rise by factors 
of 2-3 for clouds in the molecular ring of the Galaxy (Wolfire et al 2003). As discussed above 
(§2.2), depending on the galactocentric radius i?, they have densities 30 cm~^ < n < 100 cm~^, 
temperatures T ~ 50 — 100 K, and incident FUV fluxes characterized by 1< % < 4. The ratio x/'^2, 
critical to the photochemistry, varies typically from about 3-6, although we explore a somewhat 
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larger range 3-10 here. 

In the low density regime, n < 10^ cm~^, appropriate for diffuse and translucent clouds, 
there is a scaling of the results of thermochemical models. The parameters which control the 
thermochcmical structure of a cloud are just x/"-) Ccrp/n, and Ayt- We note that this scaling docs 
not apply to the denser PDR models discussed in the previous section. At high density, n > 10'^ 
cm~'^, the cooling of the gas is suppressed by collisional de-excitation of the fine structure states 
of C+ and O. Therefore, higher density gas will typically give higher gas temperatures, even when 
x/n is held constant. However, given that this scaling applies for diffuse and translucent clouds, 
we generally in this section plot our results as functions of x/i^i Ccrp/n, and Ayt- We note that 
this implies that the observations may give us a measure of Ccrp/^j but to get an estimate of Ccr-p 
we will then have to estimate the gas density n in the region observed. 



3.2.1. The total column of OH^ through a diffuse cloud of size Ayt 

Figure 10 plots the column of OH^, A^(OH+), as a function of Ayt for four values of C,crp/n2 
(recall, n2 = n/100 cm~^) and for x/"2 = 3.16. The dotted lines plot Xc(H2), the abundance of 
H2 at cloud center, and these values appear on the right of the figure. Recall Ayt is the total 
Ay through the diffuse slab, measured perpendicular to the slab and the ion columns plotted are 
perpendicular to the slab. 

Figure 10 shows that A/"(OH+) increases with Qrp/n and with Ayt, although one sees a sat- 
uration of the column with Ayt at Ayt > 0.3 for the lower three cases of Qcrpln. This saturation 
is due to the peaking of the local OH"*" abundance at lower values of Ay. Figure 10 shows that 
A^(OH+) correlates with Xc(H2), and that moderate molecular hydrogen abundances ,Tc(H2) ^0.03 
are required to obtain substantial columns of OH"'". Observed OH"*" columns in broad (Af ~ 20 km 
s~^, see §4.1.1) velocity components in the ISM are of order 3x lO-*^^ — 3x 10^^ cm~^. With this value 
of x/n2 = 3.16, we see that it requires very high cosmic ray ionization rates, Ccrp/'^2 ^3 x 10~^^ 

and high Ayt > 0.3 for a single cloud, seen face on, to produce these columns. Therefore, it 
appears that geometrical eff'ects such as seeing several clouds along the line of sight, or viewing 
the cloud obliquely, may be required to match observation. In fact, as we discuss in §4.1.1, many 
clouds may to required to cover the broad [Av ~ 20 km s~^) absorption components since single 
clouds would have much narrower velocity widths. Because of geometrical effects and variation in 
Ayt and x/i^^ the observation of the column of OH"*" is not sufficient to tightly restrict Qrp/n. 

Figure 11 plots iV(OH+) as a function of Ayt for four values of C,crp/n2 and for x/'^2 = 10: in 
other words, the same as Figure 10 but with a ratio of x/i^ that is 3.16 times higher. The main 
effect of raising x/n is to push the H/H2 transition deeper into the cloud to higher Ay. This is 
seen in the plot of Xc(H2) which rises to 0.01 at Ayt ~ 0.3 in this case, compared to Ayt '-^ 0.1 
in the previous case of x/'^2 = 3.16. Since H2 is required to make 0H+, this moves the peak of 
OH"*" to higher values of Ay, and thus pushes the rise of A^(OH"'") to higher values of Ayt than in 
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Figure 10. However, this first peak has more column of OH"*" compared with the case xl'f^2 = 3.16. 
Therefore, to produce the observed columns of 0H+ now requires only Ccrp/'n2 ^3 x 10^^^ s^^ and 
high Avt > 1 for a single cloud, seen face on. 

3.2.2. N(OR+ )/N(K20+ ) and N )/N(B.) through a diffuse cloud of size Ayt 

The results on the total column of OH"'" suggests that other observations must be brought into 
play to further restrict the cosmic ray ionization rate. Since H2O"'" is often observed as well, we 
first consider the ratio iV(0H+)/iV(H20"^). As seen in Figures 2, 3, and 5, the H20^ abundance 
peaks at higher values of Ay than the 0H+ abundance. Therefore, the ratio of the columns of 
these two ions may give us some measure of Ayt- 

Figure 12 plots the ratio A^(0H"'")/A^(H20"'") as a function of Ayt for four values of Ccrp/'^2 
and for x/''^2 = 3.16. As expected, at low values of Ayt^^-l, the ratio A^(0H"'")/A^(H20"'") is 
high, ~ 10 — 30, as we have not yet reached high enough values of Ay to include the H2O"*" peak. 
However, as Ayt increases, the ratio of the columns continues to drop until at Ayt ^ 1, the ratio is 
approximately 2-3, except in the case of the highest cosmic ray rate Ccrpln'2 = 3.16 x 10~^^ g-i|io| 
Once we have incorporated both the OH"*" and the H2O"'" abundance peaks, the column ratios are 
near unity, as noted in §3.1. 

Figure 13 plots the ratio A^(0H+)/A^(H20"'") as a function of Ayt for four values of Ccrp/'i-2 and 
for x/f^2 = 10: in other words, the same cases as Figure 12 except the ratio x/^ is raised by a factor 
of 3.16. The rise in x/i^' pushes the H/H2 transition and the peaks of OH"*" and H20^ to higher 
values of Ay. Therefore, compared to Figure 12, we see the drop in the A^(0H"'")/A^(H20+) ratio at 
higher values of Ayt. The observed ratios of A^(0H+)/A'^(H20"'") range from 3-15. If x/f^2 = 3.16 
and if we assume that Ccrp/n2 ^3 x 10~^^ s~^. Figure 12 suggests that Ayt ~ 0.1 — 0.3 to obtain 
these ratios. However, if this is true. Figure 10 implies that even with as high a value of C,crpl'f^2 = 
3 x 10^^^ s^^, we will need a "geometrical factor" of ~ 10 in order to obtain A^(OH"'")~ 10^^ cm^^. 
This geometrical factor is the combination of many clouds along the line of sight, along with the 
enhancement in column due to viewing angle of the cloud. The situation changes with x/'^2 = 10, 
as seen in Figure 13. Here, again assuming that Ccrp/n'2 ^ 3 x 10~^^ s~^, the observed ratios can be 
obtained with Ayt ~ 1 — 3. In this case, using Figure 11, the geometrical factor needs to be only 
~ 1 — 3 to produce OH"*" columns of 10^^ cm~^ when Ccrp/n'2 ~ 3 x 10^^^ s^^. 

Observations of HI 21 cm are often available along the same lines of sight as the OH"^ and 



^''The ("crp/n2 = 3.16 x 10~^^ case is extreme; cosmic ray ionization of H produces an electron abundance of 
Xe — 3x 10~^ . In this case, electrons and not H2 dominate the destruction of OH^ and H20^ even when a::(H2) ~ 0.5. 
Therefore, the fraction of cosmic ray ionizations that produce OH"*" is higher than the fraction that produce H20^, 
even at the peaks of their abundance. Hence we have much more OH"'" than H20"''as seen in this figure, and, once 
we reach high enough Ayt where the ion chemistry shown in Figure l(top) dominates, the ratio rises as a function of 
Avt- This holds as well in X-ray ionized regions with high ionization rates. 
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H2O''' measurements. In this case, the column of HI, A'^(II), along the line of sight associated 
with the velocity feature of the ions can be estimated. To be precise, the observations directly give 
A^(H)/T, where T is an average temperature along the line of sight. In Figures 14 and 15 we plot the 
ratio iV(0H+)/[A^(H)/r2] on the vertical axis {T2 = T/lOO K) and the ratio A^(0H+)/Ar(H20+) 
on the horizontal axis for two fixed values of = 3.16 and 10. The results of our 2 sided diffuse 
cloud models are shown as contours on this figure. As noted above, higher average abundances of 
0H+ (or iV(OH+)/iV(H)) requires higher Ccrp/n. On the other hand, higher 7V(OH+)/7V(H20+) 
requires either low Ayt or high x/^- the next section we apply Figures 10, 11, 14, and 15 to the 
observational data to constrain Ccrp/n' along sightlines to W49N. 



3. 2. 3. The possible effects of turbulence 

We have run additional models to estimate how turbulence might affect the production of 
0H+, 1120"*", and HsO^. It has been recognized for several decades that standard astrochemical 
models for diffuse molecular clouds greatly underpredict the observed column densities (~ 10^^ 
cm~^) of CH"*" (e.g., Falgarone et al 2010a, b; and references therein.) The models presented here 
are no exception; in our models with x ~ 1, n ~ 100 cm~^, and Ayt > 1, for example, we obtain 
a predicted CH"*" column density ~ 3 x lO^'' cm~^, more than two orders of magnitude lower than 
that typically observed. Both shock waves with large scale order (Elitzur & Watson 1978, Draine & 
Katz 1986, Flower &: Pineau des Forets 1998), and interstellar turbulence (e.g. Joulain et al. 1998; 
Lesaffre et al 2007, Godard et al 2009, Falgarone et al 2010a,b) have been proposed as sources of 
heating that could enhance the CH+ abundance within some small fraction of the cloud volume; 
in these models, CH+ is formed by the endothermic reaction of C"^ with H2. It has also been 
long recognized (e.g. Draine &: Katz 1986) that models invoking ion-neutral drift - either in C-type 
shock waves or in turbulent dissipation regions - are most successful in simultaneously matching the 
observed column densities of CH"*" and OH. In particular, while models in which the gas temperature 
is merely elevated without significant ion- neutral drift tend to overpredict the OH/CH"*" ratio (the 
neutral-neutral endothermic reaction O + H2 — > OH + H being enhanced along with the reaction 
of C"^ with H2), models incorporating ion- neutral drift preferentially enhance the endothermic 
ion- molecule reactions that produce CH"*" (and SH+; Godard et al. 2012), without producing 
more OH than the observations permit. We have crudely estimated the effects of turbulence by 
positing an enhanced rate of endothermic ion-neutral reactions in a small fraction of the cloud 
volume. Motivated by recent calculations performed by Myers & McKee (private communication), 
we adopted enhanced ion-neutral reaction rates - equal to the thermal rate coefficients at 1000 
K - in 3% of the cloud volume: this leads to a predicted CH"*" column in accord with what is 
typically observed^ We found that turbulence does not greatly enhance the OH"*", H20^, and 



^^In particular, we add 97% of the columns from our standard model to 3% of our columns from the same model 
but with the enhanced ion-neutral rates to obtain the total columns of, for example, CH"*". 
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HaO"*" column densities through such a cloud. There was a modest enhancement of these oxygen 
hydride ions within the region of ion- neutral drift, due to the faster reaction of with O, but 
since that region occupies only a small fraction of the total volume, this effect did not affect the 
total column densities significantly. Therefore, we conclude that turbulent dissipation or shocks in 
small regions of the clouds are unlikely to affect our conclusions concerning 0H+, H2O+ or H3O+ 
in the models presented here. 

4. OBSERVATIONS AND COMPARISON WITH MODEL PREDICTIONS 

4.1. OH+, H2O+, and H3O+ Absorption in Diffuse Galactic Molecular Clouds 

As discussed in §1, recent ifersc/ie/ observations of bright Galactic continuum sources have re- 
vealed OH"*" and H20'^ absorption features arising in multiple foreground clouds along the targeted 
sight- lines (Gerin et al. 2010; NIO; Neufeld et al. 2011). The three crosses in Figure 14 denote 
the A/"(0H+)/Ar(H20+) and N{0R+)/[N{R)/T2] ratios for the two such sources for which results 
have been reported to date: G10.6-0.4 (a.k.a. W31C) (the lower limits) and W49N (2 points to 
the right). Those given for G10.6-0.4 apply to all the material along the sight-line, while those 
given for W49N apply to two separate velocity intervals (30-50 km s^^ and 50-78 km s~^). The 
OH"*" absorption in GlO.6-0.4 is largely saturated, and thus the plotted values are lower limits. For 
W49N, the columns of H derived assumed T = 100 K or T2 = 1; the quoted values of N{R)/T2 
are 6.95 x 10^^ cm^^ for the 30 - 50 km s~^ feature and 7.23 x 10^^ cm^^ for the 50 - 78 km s~^ 
feature. The total columns A''(OH^) observed in these two velocity intervals are 3.6 x 10^^ and 
2.1 x 10^^ cm~^, respectively. 

The A^(0H+)/A^(H20'^) ratios measured toward W49N (~ 10) and the wide velocity range 
of the absorbing clouds suggest that the absorbing material is comprised of multiple clouds of 
small extinction. The line of sight to W49N is approximately 11 kpc long, and passes as close as 
R = 5 kpc to the Galactic center. Assuming x/"2 = 3.16, models with Ayt = 0.32 and 0.25 best 
account for all the data (Figure 14), with cosmic-ray ionization rates of ~ 6 and 4 x 10~^^n2S~^, 
respectively, for the 30 — 50 km s~^ and 50 — 78 km s~^ velocity intervals. As seen in Figure 10 for 
x/n'2 = 3.16, the column of 0H+ in a single cloud of size Ayt = 0.32 is A'"(OH+)~ 2 x 10^^ cm~^ 
and for Ayt = 0.25 is A^(OH+)~ 1.5 x 10"*^'^ cm^^. Therefore, we either need 15 clouds along 
the line of sight, or a smaller number but with geometrical enhancement effects. Note that a single 
cloud might have an OH"*" absorption linewidth of only ~ 1 — 3 km s~^. Therefore, to cover the 
broad absorption features we are modeling {Av = 20 or 28 km s~^), we need 10-30 clouds spread 
out along the line of sight so that their galactic rotational velocities coupled with their individual 
line widths cover this velocity range. Therefore, this x/"^2 = 3.16 model has barely enough clouds 
to produce the relatively smooth absorption feature observed. Assuming x/''^2 = 10, models with 
Avt — 0.75 and Ayt — 0.63 and with Ccrp/?^2 ~ 2 x 10~^^ and 1.2 x 10^^^ s^^ respectively, best 
account for the data (see Figure 15). With these values of Ayt and Ccrp/^2, inspection of Figure 
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11 shows that we require ~ 12 clouds along the line of sight. This value of has even fewer 
clouds to cover the broad absorption hue observed. Moreover, the temperatures in clouds with 

= 10 are higher (~ 200 K) than is typically observed near the solar neighborhood (Heiles &: 
Troland 2003). And finally, the Wolfire et al (2003) results suggest that along the line of sight to 
W49N, the average ~ 4. A similar modeling of the case x/n2 = 1, not shown in the figures, 

gives C,crp/n2 = 2.5 x 10^^^ and 2.0 x 10~^^ s~^ with Ayt = 0.16 and 0.08, respectively. Again, 
because of the high cosmic ray rates, the gas temperatures are high in these clouds (~ 150 K). In 
addition, these cosmic ray rates seem improbably high. This model is driven to high cosmic ray 
rates in part because otherwise (with low the gas is cold, which reduces the rate that H"*" can 

charge exchange with O. Therefore, we finally conclude that our models suggest typical Ayt ~ 0.3, 

~ 3, and Ccrp/'T'2 ~ 6 x 10~^^ s~^ for the 30 — 50 km s^^ component, and C,crp/n2 ~ 4 x 10"^^ 
s~^ for the 50 — 78 km s~^ component. 

A better approach is to ask, "What distribution of cloud Ayt or hydrogen nucleus column N 
through the cloud will we encounter in traversing the ~ 11 kpc to W49N?" HI 21 cm observational 
data suggest that dud/dN oc A^~^, where Ud is the number of clouds for N > 2.6 x 10^'^ cm~^, and 
proportional to A^~^ for ^ 2.6 x lO^'^ cm" -2 (Heiles & Troland 2005)0 We have crudely integrated 
a distribution of clouds with this dependence, with the constant of proportionality derived by 
requiring the integral to obtain a total column of A(H)= 7 x 10^^ cm~^, the observed atomic H 
column in each velocity component of W49N. We use A = 6 x 10^^ cm~^ as an upper limit for 
the integration; above this, we enter small number statistics since the total column is of this order. 
The lower limit does not enter the integration significantly, as long as it is << 2.6 x 10^'^ cm~^, 
since these clouds contain very little of the total column of any species. The integration removes 
Ayt as a variable, and leaves only the parameters C,crp/n2 and x/^2 as free parameters to match 
the observed 0H"''/H20^ ratio and the total column of OH"*". We find that with x/^2 = 3.16 and 
Ccrpln2 = 3.7 X 10"^^ s-\ we obtain A(OH+)~ 2.1 x 10^"^ cm^^ and Ar(H20+)~ 2.4 x 10^^ cm-^. 
These values match the velocity interval 50 — 78 km s~^ extremely well. Increasing the cosmic ray 
rate to Ccrp/A'2 = 6.3 x 10"^^ s-\ we obtain A(OH+)~ 3.6 x 10^^ cm'^ and A(H20+)~ 4.0 x lO^^ 
cm~-^, a very good fit to the lower velocity component. Considering the small number statistics, 
especially of the higher Ayt clouds, this is very good agreement. Interestingly, the main contribution 
to A(OH+) comes from near Ayt = 0.2 — 0.6. This is because of the dependence of A(OH+) on 
Ayt (see Figure 10). Note that A(OH+) rises sharply with Ayt until about Ayt = 0.3, and then 
it levels off. Therefore, the integral is weighted to that size cloud in the range of the turnover 
{Ayt ~ 0.3 — 0.6), which has the peak OH"*" abundance (averaged through the cloud). In addition, 
the cloud distribution steepens when Ayt > 0.13, thereby decreasing the contribution from higher 
Ayt clouds. The total number of clouds in the range Ayt ~ 0.1 — 1 is approximately 20. We note 
that the H20"'" columns have greater contribution from higher Ayt clouds, which are less numerous, 
and therefore we predict the H20^ absorption feature to have greater fluctuations. 



Recall that A'^ is the column of H nuclei; however, in the above range, generally ~ A^(H). 
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We next examine the cases xl'^i = 1 ^^i^fl 10- Assuming xl'^i — 10) we find that we need 
cosmic ray rates of 1.1 x 10^"*^^ s^^ and 7 x 10^^*^ s^^ to match the 0H+ columns in the two velocity 
components, respectively. However, the 0H+/H20'^ column ratio is 12.1, somewhat higher than 
observed, and again, the clouds are too warm (~ 200 K). The high field pushes the peaks of 0H+ 
and H20'^ deeper into the cloud, and our clouds tend to truncate the H2O''" column, leading to 
the high OH+/H2O'*" ratio. We need much higher cosmic ray rates than our solution above for 
— 10 with clouds of single Ayi ~ 0.7 because the cloud distribution leads to a significant 
column of H from clouds of low Ayi which have very little 0H+. Assuming x/'^2 = 1) wc find 
that we need cosmic ray rates of 2 x 10~^^ s~^ and 1.2 x 10~^^ s~^ to match the OH^ columns 
in the two velocity components, respectively. However, the OH+/H2O+ column ratio is 5.6, lower 
than observed. With low x/n2, the clouds tend to contain both 0H+ and H2O"'" peaks, driving 
the ratio down. We conclude that the best fit is with x/n2 = 3.16, giving C,crp/n2 ~ 6 x 10~^^ 
s~^ and 4 x 10~^^ s~^ for the two velocity components, in good agreement with our analysis above 
that did not use the cloud distribution. However, the cloud distribution gives added weight to the 
conclusion that x/ra2 ~ 3. 

The cosmic ray rates derived from our models are roughly 2-4 times larger than those inferred 
by NIO, who used a simple analytic treatment to infer cosmic ray ionization rates of 0.6 and 
1.2 x 10~^^ (n2/e)s~^, where e is the fraction of cosmic-ray ionizations that lead to 0H+. Using 
pure gas-phase model results from the Meudon PDR code (Le Petit et al. 2006), NIO found that 
e lay in the range 0.5 - 1.0 for a wide variety of cloud conditions; this would imply cosmic ray 
ionization rates in the range ~ 0.6 — 2.4 x 10~^^ n2 s~^ for the material along the sight-line. The 
main source of the difference between our models and those of NIO is our inclusion of PAH^ - H"*" 
recombination and the charge exchange of H+ with neutral PAH in the present work, processes - 
previously omitted - that reduce the fraction of cosmic-ray ionizations that lead to 0H+ production 
(in our models e ~ 0.1 — 0.3 in the regions where most of the 0H+ lies). 

Figures 16 and 17 repeat Figures 14 and 15 except that the rates of PAH and PAH~ neu- 
tralization of H"*" have been reduced by > 4, so that they are negligible compared to the radiative 
recombination of H+ with electrons (see Appendix B) . Note that this reduction could be achieved 
by either lower PAH abundances or lower PAH rate coefficients with H^, or a combination of both. 
Another possibility is increased photoionization rates of PAH^ and PAH. Figure 16 shows that the 
best fit to the data for xjn = 3.16 is Ayi ^ 0.316 and C,crp/n2 1 — 2 x 10~^^ s~^, in complete 
agreement with NIO. Therefore, the PAH chemistry is extremely important in deriving cosmic ray 
rates from the 0H+ and H2O+ ions. We consider the rates from Figures 16 and 17 lower limits to 
the cosmic ray rates, whereas the rates from Figures 14 and 15 could be considered upper limits, 
although they represent our best estimates of PAH chemistry at this time. 

The cosmic-ray ionization rates suggested by the predictions shown in Figure 14 are also a 
factor ~ 5 — 10 larger than those typically inferred from observations of Hg" in the Galactic disk 
(Indriolo et al. 2007, Indriolo & McCall 2012). Indriolo & McCall conclude, however, that cosmic 
ray rates vary by more than an order of magnitude over various sight lines, from Qcrp^ 0.7 x 10~^^ 
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to 4.6 X 10~^^ s~^, with a mean of ~ 1.5 x lO"-*^^ (note that we have converted their 
total ionization rate per H2 to our primary ionization rate per H nucleus by dividing their values 
by 2.3; we also note that their typical density was n ^ 200 cm^^ so their results imply Ccrp/'^2 
that are two times smaller than these values. Their highest values of Ccrp/n2 are still a factor of 
~ 2 below what we find for W49N, so there appears to be a small discrepancy. We point out 
that the cosmic ray rate derived from H3" observations is directly proportional to the assumed 
electron abundance. Indriolo & McCall assumed that the electrons were supplied by C"*", with an 
abundance of 1.5 x 10"*^. However, in our models of diffuse clouds, where typically the gas density 
is ~ 100 cm~^, the cosmic ray ionization of H can produce comparable or greater abundances of 
H"*" than C"*". Therefore, we find for Ccrp/n2 = 2 x 10~^^ s~^, for example, that — 3 x 10~^ when 
a;(H2)~ 0.1 and Xe — 2 x 10~^ when x(H2)~ 0.4 . This implies that the Indriolo et al rates may 
need to be increased, by factors as much as ~ 2, depending on Ccrp/?^2- Given the small number 
of sources toward which A'^(0H"'")/iV(H20+) and iV(OH^)/A^(H) ratios have so far been reported, 
the proven variation in cosmic ray rates along different lines of sight, the lack of a common sight 
line to compare the cosmic ray rates derived by versus by OH+and H2O+, the possibility of 
higher electron abundances than assumed by Indriolo and coworkers, and the lack of certainty in 
PAH chemistry, it is not yet possible to draw firm conclusions about whether our results differ 
significantly from that of Indriolo et al. What is clear is that cosmic ray rates in the diffuse ISM 
are higher than was previously thought. 

Herschel has also detected absorption by H3O+ in foreground material along the sight-line to 
G10.6 - 0.4 (Gerin et al. 2010). The average H3O+/H2O+ ratio along the sight-line is 0.7, but the 
distribution of HaO^ is clearly different from that of H20^ and OH^, with the HsO^ concentrated 
within a single narrow velocity component where the H2 fraction is presumably largest. This cloud 
is unassociated with the source, but is probably the foreground cloud with the largest Ay along 
the line of sight. It appears in the HF spectrum (Neufeld et al. 2010) and CO spectrum (C. Vastel, 
private communication; see also Corbel & Eikenberry 2004), suggesting a high molecular fraction 
and therefore substantial Ay- The HaO"*" column in this narrow component is ~ 1 — 2 x 10^^ cm~^, 
consistent with our predictions for the second peak in a high Ay cloud. 



4.2. H3O"'" Emission from Dense Galactic Molecular Clouds 

PvDK92 searched for the 396, 364, and 307 GHz inversion-rotation lines of H3O+ in 20 Galactic 
targets - including GMCs, star forming regions. Galactic center sources and a few evolved stars - 
and reported definitive detections of H3O+ emission from the G34. 3-1-0. 15 hot core, two positions 
in Sgr B2 (a Galactic center GMC), and the W3 IRS5, W3(0H), W51M and Orion KL high- 
mass star forming regions. In the case of Orion and Sgr B2, H3O+ emissions had been previously 
discovered by Wootten et al. (1991). The HsO"*" column densities derived by PvDK92 for the high- 
mass cores W3 1RS5, W3(0H), W51M are fairly similar, ranging from a few xlO^"^ to a few xlO^"^ 
cm~^, and arguing for a similar environment and formation mechanism. Given the large critical 
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densities for the observed transitions, PvDK92 inferred relatively high densities (n > 10^ cm~^) 
and temperatures (T > 50 K) for the HsO"*" emitting gas. The HsO"*" column densities observed 
in these high-mass cores are in good agreement with the predictions of our models for either the 
case in which PAHs and/or very small grains are assumed to be present or the case where PAHs 
and low ionization potential metals are highly depletecQ); the PAH case predicts HsO"*" column 
densities of 3 x 10^^ - lO^"^ cm ^, regardless of density n or C,crp- 

Herschel's HIFI instrument has allowed the detection of far-infrared and submillimeter HsO^ 
lines that are not accessible to ground-based observatories. Several H3O+ emission lines in the 
THz domain have been reported towards high-mass YSOs in the W3 IRS5 region (Benz et al. 
2010). A rotation diagram combining HIFI data with lower frequency data from PvDK92 indicates 
a rotational temperature T^ot ~ 239 K, suggesting that the observed emission may arise either from 
hot and dense gas or is radiatively excited by continuum photons from hot dust grains. Benz et 
al. (2010) inferred A^(H30"^) = 8.5 it 2 x 10^'^ cm~^ and interpreted the HsO"*" emission as arising 
from the outflow walls heated and irradiated by the FUV radiation field from massive YSOs. Their 
inferred HsO^ column densities and rotational temperatures are consistent with our models of dense 
gas illuminated by a relatively strong FUV field (where large kinetic temperatures are achieved at 
relatively low Ay ~ 1 — 2). 



4.3. Combined H3O+ Emission and Absorption in Strong Far-infrared Continuum 

Sources 

For sources where a strong infrared radiation field is present, submillimeter HaO"^ line emission 
can be accompanied by far-infrared absorption within an extended envelope. The H3O''' Jk='^i — 
and 2q — 1q inversion-rotation and lj~ — 1^ pure-inversion FIR absorption lines at 100.577, 100.869 
and 181.054 /um were detected by ISO towards the optically thick FIR continuum emission of Sgr 
B2 (Goicoechea &: Cernicharo 2001). An H3O"'' column density of ~ 1.6 x 10^"^ cm~^ was inferred 
in the warm and low density extended envelope of Sgr B2 (n = 10'^ — 10^ cm~^, x = 10^ ~ 10^ 
from [01] and [CII] observations; Goicoechea et al. 2004). However, the p-HsO"^ 3^ — 22 364 
GHz emission line was subsequently mapped around the Sgr B2(M) core at 18" angular resolution 
with APEX by van der Tak et al. (2006), who inferred column densities ~ 10^^ — 10^^ cm~^ on 
the basis of a detailed excitation model. This value is much larger than the estimates obtained 
from FIR absorption lines, presumably because the latter only probe the outer envelope of the 
source whereas the submillimeter observations probe material to a much larger depth. The HaO"*" 
column densities inferred by van der Tak et al. (2006) are much larger than our model predictions 



^^If PAHs and very small grains coagulate into larger aggregates deeper inside the cloud (e.g., Boulanger et al. 
1990; Rapacioli et al. 2006) or onto larger grain surfaces, our models would predict H3O''' column densities that 
are roughly an order of magnitude lower than the values inferred from the above observations, unless low ionization 
potential metal atoms are highly depleted at the second peak of HsO"*". 
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for externally-illuminated clouds; this discrepancy likely reflects the presence of a strong internal 
luminosity source in Sgr B2 that heats the dust grains within the interior and prevents the freeze-out 
of oxygen nuclei (see Appendix C). 

4.4. Extragalactic, OH+, H2O+ and H3O+ 

Thanks to the much improved sensitivity of space- and ground-based receivers and detectors, 
molecular ions can now be detected outside the Milky Way. Extragalactic HaO"*" was first detected 
through JCMT observations of the P-H3O+ 364 GHz emission line towards the prototypical ul- 
traluminous infrared galaxy Arp 220, and towards M82, an evolved starburst (van der Tak et al. 
2008). H3O+ was subsequently detected towards IC342, NGC253, NGC1068, NGC4418 and NGC 
6240, and upper limits obtained towards IRAS 15250 and Arp 299 galaxies (Aalto et al. 2011). 
Except for IC342 and M82, the typical H3O+ column densities (~ 10^^ - 10^'' cm ) derived from 
extragalactic H3O"'" detections are much larger than the predictions of our models for single clouds. 
These observations are of HsO"*" in emission and the authors assume T^x ~ T^ad to obtain these 
columns with error bars of factor of 2. As in the case of Sgr B2, the discrepancy between our 
models and the observations may reflect the effect of enhanced dust temperatures that prevent the 
freeze-out of elemental oxygen. Alternatively, there may be a large number of clouds along the line 
of sight. 

For the starburst galaxies M82 and IC342 , however, where H30^ column densities ('^ 10"^^ 
cm^^) are inferred from the observations, enhanced dust temperatures or ionization rates are not 
apparently required. For example, van der Tak et al (2008) inferred Af(H30''~)~ 1.1 x 10^^ cm~^ 
in M82 from an emission line which contained a combination of a broad (Aw ~ 260 km s~^) and a 
narrow (Av ~ 40 km s~^) component. In addition, a recent Herschel/llWl detection (Weiss et al. 
2010) of the 0-II2O"'" 111 — Ooo ground-state line in absorption towards M82 revealed A'"(H20"'")~ 
2.9 X 10"*^^ cm^^, arising from a Av ~ 77 km s^^ feature. The large width of the lines suggest 
numerous clouds in the beam. Overall, these numbers are consistent with our models with the 
1120"*" arising from the flrst (and possibly second) peak of numerous relatively low density clouds, 
whereas the H3O+ arises from the second peaks of many molecular clouds in the beam. 

5. Summary and Conclusions 

We model the production of 0H+, 112©"^, and HsO"^ in interstellar clouds, using a steady state 
photodissociation region code that treats the freeze-out of gas species, grain surface chemistry, and 
desorption of ices from grains. The code includes PAHs, which have an important effect on the 
chemistry. 

As a function of depth or Ay into a molecular cloud, the ions tend to have two peaks, one 
at low Ay < 1 and one at high Ay ~ 6, the exact value depending on x/n. In the first peak the 
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ions are created by the cosmic ray ionization of H to H"*" , followed by reactions with O and H2 (see 
Figure 1 top). These peaks appear on the surfaces of molecular clouds as well as in diffuse clouds. 
PAHs can lower the production of the ions here, by neutralizing and interrupting the reaction 
chain (Figure 1 top). At most, they depress the ion abundances and columns by a factor of 3. 

In molecular clouds a significant fraction of the column density of OH+and H20^ is found in 
the first peak at the surface {Ay < 1) of the cloud. For relatively low values of the incident far 
ultraviolet flux on the cloud (%< 1000), the columns of OH"'" and H2O+ scale as the cosmic ray 
ionization rate divided by the gas density. Roughly, the columns of 0H+ , H2O+, and H3O+ in the 
first peak arc 2.2 x W^'^x^^^iCcrp/'^ x lO^^^ s^^) (100 cm-^/n) cm'^, 1.5 X lO^^x^-^^ (Ccrp/2 X 10""^^ 
s"^) (100 cm-^/n) cm'^, and 9 x 10^^ (Ccrp/2 x 10"^^ s'^) (100 cm-^/n) cm'^. 

There is a second peak in 0H+, H2O''", and H3O+ abundances at larger depths {Ay ^ 6) in 
molecular clouds, when the second route to OH"'" formation, initiated by the cosmic ray ionization of 
H2 becomes dominant (Figure 1 bottom). Here, lower electron abundances enhance the abundances 
of the ions by lowering the electronic recombination rate of H^ (which raises the abundance of H3" 
and the formation rates of the three ions), and by lowering the electronic recombination of H3O''" 
(its main destruction path). However, even deeper in the cloud, the oxygen freezes out as water 
ice on the grain surfaces, and the ion abundances fall as their formation rates fall, being starved 
for gas phase elemental oxygen. If PAHs or VSGs are present at these high values of Ay ^ 6, 
the electron abundance at the second peak is controlled by association with neutral PAH or VSG. 
In this case, rather surprisingly, the column of H3O"'" (~ 4 x 10^'^ cm~'^) in the second peak and 
the peak abundance (~ 10""^) of H3O"'" is independent of both Ccrp and n. Raising the cosmic ray 
ionization rate increases the production rate of H3O"'", but it also increases the destruction rate by 
electrons by the same amount. If PAHs and VSGs are not present at high Ay, the column of H3O"'" 
depends mainly on the gas phase abundance of elemental S, Si, Mg. and Fe at the second peak. 
If these abundances are low, < 3 x 10~^, then columns comparable to the PAH case arc obtained. 
However, for higher abundances of the metals, the H3O"'" column is reduced by a factor of roughly 
10. The columns of OH"*" and H2O"'" in the second peaks are usually of order of the columns in the 
first peaks. 

The observed H3O"'" columns of ^ 4 x 10^^ cm~^ seen in many dense Galactic molecular clouds 
therefore imply that either PAHs or VSGs are present deep in molecular clouds or that the depletion 
of PAHs and VSGs are accompanied by a very significant depletion of gas phase metals. There are 
a few observations that imply much higher column (~ 10^^ — 10^^ cm~^) of H3O"'", and these can 
only be explained in the context of our models as arising in very high Ay clouds where high gas 
phase elemental O abundances are present throughout due to either desorption processes or time 
dependent effects. We suspect that the grains in these sources may be so warm, ^ 100 K, that 
thermal desorption keeps water ice from depleting the oxygen reservoir. In the case of extragalactic 
sources with very broad velocity features, the columns may be increased by the presence of a large 
number of clouds along the line of sight. 
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For high values of the incident far ultraviolet flux {x ^ 1000) and high gas densities ( ^ 10'* 
cm~^), producing warm (T > 300 K) gas with significant H2 abundance at Ay ~ 1 — 2, chemical 
reactions initiated by the photoionization of carbon in the mainly atomic surface can form ionized 
hydrogen, which then leads to the formation of 0H+, H2O+, and HsO"^. In this case, their columns 
(typically, of order 3 x 10^^ cm~-^) are not related to the cosmic ray ionization rates. HsO"*" emission 
from W3 IRS 5 may be an example of such conditions. 

We also model diffuse and translucent clouds in the interstellar medium, and show how obser- 
vations of iV(OH+)/iV(H) (typically lO"*^ - lO""^) and 7V(OH+)/iV(H20+) (typically 3-15) can be 
used to estimate Ccrp/n. The ratio A^(OH^)/A^(H), which is essentially the average abundance of 
OH^ in all the clouds along the line of sight (within the same absorption velocity component), is 
mainly a measure of Ccrp/n. The ratio of the OH"*" column to the 112©"^ column in diffuse clouds 
is mainly dependent on the ratio x/"- and Ayt (i-e, the total hydrogen column through a single 
typical cloud). If x/'t- is known, or at least constrained to a narrow range such as x/'^2 ~ 3 — 6 
typical of diffuse clouds in random locations along the line of sight in the Galaxy, then observations 
of iV(OH+)/iV(H), iV(0H+)/iV(H20+), and iV(OH+) can determine Ccrp/n, Ayt , and a geomet- 
rical factor which is a combination of the number of clouds along the line of sight and the typical 
angle that the line of sight makes as it passes through each diffuse cloud slab. Using the observed 
distribution of Ayt in diffuse clouds, the models can provide an estimate of x/n2- 

We discuss the relation of our models to recent observations of OH"*" and 1120"*" by the Herschel 
Space Observatory, and the ability to constrain the cosmic ray ionization rate through comparison 
of observations with these models. We conclude that the cosmic ray primary ionization rates 
Ccrp in the observed foreground diffuse clouds towards W49N have values of approximately Ccrp~ 
4 — 6 X 10^^®(n/100 cm^'^) s^^, if our adopted PAH chemistry is correct0 We find a hard lower 
limit of Ccrp/n2 ^ 1 — 2 x 10^"*^^ s^^ by neglecting PAH chemistry in W49N. Our best fit models 
suggest that x/'^2 ~ 3 in the diffuse clouds towards W49N. Our W49 models suggest that, in terms 
of producing 0H+ and H2O"'" column, the typical Ayt ~ 0.3 through a single cloud towards W49N. 
and that a diffuse cloud Ayt distribution measured by Heiles & Troland (2005) fits the data very 
well. To produce the total column of OH"*" and H observed in the two velocity components requires 
~ 20 clouds in each component along the line of sight. We discuss differences between our estimates 
and those of NIO and Indriolo et al (2007), pointing out the former neglected PAH chemistry and 
the latter may have slightly underestimated the electron abundance in the diffuse foreground clouds 
with the highest cosmic ray ionization rates. We look forward to further observations of OH"*" and 
H2O"'" along many sight lines to probe the cosmic ray ionization rates throughout the Galaxy. 
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Appendix A: Tables of rate coefficients and adopted abundances 

Table 1 presents some of the key reactions and/or some of the reaction rates that have changed 
since H09. We emphasize that this is not a complete listing of the ~ 300 chemical reactions in 

the PDR code described in HOD. The PAH rates arc scaled by the factor $paH) a scaling factor 
introduced in Wolfirc ct al (2008) that incorporates the uncertainties in PAH reaction rates, PAH 
sizes, and PAH abundances. Wolfire et al (2008) found from a semi-empirical fit of our PDR models 
to C, C+, H, and H2 column densities in diffuse clouds that $pah ~ 0.5, and we have adopted that 
value in all of our PDR models. 

Table 2 presents the standard gas phase abundances and grain properties that we have adopted 
in most of our PDR models. As discussed in the text, we have also run models with these values 
changed to test the sensitivity of the results to the standard values. In particular, we have run 
models with x(PAH)= 0; with x{Si), x{Fe), x{S), and x(Mg) all set to 10""^ or 10~^; and with x(0) 
= 4.5 X 10-^. 
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Table 1. Reaction Rates 



Reaction 



Rate Coefficient 



PAH- + H+ 
PAH° + H+ 
PAH+ + e - 
PARO + e ^ 
PAR- + C+ 

paro + C+ 
PAR° + hu - 
PAR- + hv 
PARO + hu - 
PAR- + hv 
C + hu^C 
R2O + hu^, 
R2O + hi^^, 
R2O + hi^^, 
OR + /iz/ - 
OR + hu- 
0R+ + hiy 
CR+ + hu 
R + CR^ 
R2 + CR - 
R+ + e ^ 
R^ + e ^ 
Rj + e ^ 
R^ + e ^ 
0R+ + e - 
R2O+ + e 
R2O+ + e 
R2O+ + e 
R3O+ + e 
R3O+ + e 
R3O+ + e 
R3O+ + e 
R+ + CO 



PARO + R 
PAR+ + R 
PARO 
PAR- 

^ PARO + C 
^ PAR+ + C 
PAR+ + e 
PAR° + e 
PAR+ + e 
^ PARO + e 

> R2O+ + e 
^ OR + R 
^0 + R2 
+ R 
OR+ + e 
^ 0+ + R 
^R+ + C 
■R+ + e 
->R^ + e 
R 

R + R 
R2 + R 
R + H + R 

+ R 
^R2 + 
^ OR + R 
^ + R + R 
^ OR + R + R 
^ R2O + R 
^ R2 + OR 
^ R2 + O + R 
RCO+ + R2 



8.1 X 10 
7.0 X 10 

3.4 X 10 

3.0 X 10 
2.3 X 10 

2.0 X 10 
2.8 X 10 
5.7 X 10 

3.5 X 10 

1.7 X 10 

3.1 X 10 

3.1 X 10 
7.5 X 10 

4.8 X 10 

3.9 X 10 

2.2 X 10 
1.1 X 10 

3.3 X 10 



-^$PAH(r/300 K)-0-50 cm3 s- 
-^$PAH cm^ s 



-1 a 



■5^>pah(T/300 K)-0-™ cm3 



-6 



-7 



-1 a 



$PAH(r/300 K) 



-0.50 



1 a 



-1 a 



-1 a,b 



^$PAH cm^ s ^ ^''^ 




^Xexp(- 


-2.34Av) s- 


1 c,d 


^Xexp(- 


-1.09Av) s- 


1 c,e 


®Xexp(- 


-2.45^v) s- 


1 f,d 


^Xexp(- 


-1.77^v) s- 


1 f,e 


i°Xexp( 


-3.33Av) s 


-1 g 


i^Xexpl 


-3.90Av) s 


-1 g 


exp( 


-1.70^v) s 


-1 S,h 


^^Xexp( 


-2.20^v) s 


-1 S,h 


^°Xexp( 


-1.70Av) s 


-1 g,h 


^^Xexp( 


-4.05Av) s 


-1 i 


^^Xexp( 


-3.50Av) s 


-1 g 


^°Xexp( 


-2.94Av) s 


-1 g 



c 

2C 



crp 



crp 



-1 j 
j-1 j 



3.5 X 10" 


-i2(T/300 K) 


0.75 




-1 k 


1.6 X 10" 


-^(r/300 K)- 


0.43 






-1 k,l 


3.4 X 10" 


-^(r/300 K)- 


0.50 


cm^ 


3- 


-1 m 


3.4 X 10' 


-^(r/300 K)- 


0.50 


cm^ 


s" 


-1 m 


3.8 X 10- 


-^(r/300 K)- 


-0.50 


cm^ 


s" 


-1 k,l 


3.9 X 10- 


-^(r/300 K)- 


0.50 


3 

cm 




-1 k 


8.6 X 10- 


"S(r/300 K)- 


0.50 


3 

cm 




-1 k 


3.1 X 10- 


-^(r/300 K)- 


0.50 


3 

cm 


g- 


-1 k 


3.4 X 10- 


-^(r/300 K)- 


0.74 


cm^ 


3- 


-1 n 


1.4 X 10- 


"^(r/300 K)- 


0.74 


cm^ 


s' 


-1 n 


7.9 X 10- 


-8(r/300 K)- 


■0.74 


cm^ 


3- 


-1 n 


7.4 X 10- 


-9(r/300 K)- 


0.74 


cm^ 


s' 


-1 n 


1.7 X 10- 


-9 gj^3 g-l k 
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Table 1 — Continued 



Reaction 



Rate Coefficient 












C+ + OH^CO+ + H 


2.9 X 


10" 


9(r/300 K)-°-33 cm^ P 


C+ + H2 ^ CH+ + H 


1.0 X 


10- 


10 exp(-r/4640 K) cm^ s'^ ^ 


C0+ + H ^ CO + H+ 


7.5 X 


10- 


cm-^ 3-^ ^ 


+ H+ ^ 0+ + H 


q 






+ H ^ OH + /izy 


9.9 X 


10" 


19(r/300 K)-0-38 k 


0+ + H ^ H+ + 


5.7 X 


10" 


10(r/300 K)-0-36e8-6K/T g-1 q 


0+ + H2 ^ 0H+ + H 


1.7 X 


10- 


^ cm^ ^ 


0H+ + H2 ^ H2O+ + H 


1.0 X 


10- 




H2O+ + H2 ^ H3O+ + H 


6.4 X 


10- 


cm^ 3-^ ^ 


H+ + H2 ^ H+ + H 


2.1 X 


10- 




H+ + ^ 0H+ + H2 


8.4 X 


10" 


cm^ ^ 


H+ + ^ H2O+ + H 


3.6 X 


10" 


cm^ 3~^ ^ 


H2 + ^ OH + H 


3.40 > 


< 10 


-i3(r/300 K)2-6^exp(-r/3160 K) cm^ 



-,-1 k 



H2O"'" ha3 no photodi330ciation channels longward of 13. 61 eV and the rate is set to zero 

(jvan Dishoeck et al.ll2006l ). 



'^Non-photo PAH rates are calculated using the equations of Draiiie Sutin J 19871). R epresen- 
tative rates are given at T = 300 K for disk PAHs. $pah = 0.5 from Iwolfire et al.l (|2008l ). 

'^Additional collisonal rates scale as (m)"*^'^ where m is the mass of the collision partner. 



is the FUV field measured in units of the iDraind (|l978l ) field. Rate for Nc = 100. The shape 
of the FUV and optical field from Mathis et al. ( 19831 ) used for x = 1- 



'^Absorption cross section and ionization potentia l of circumovalene (IP=5.7 eV; iMalloci et al 



20071 ). and linear yield function ([Jochims et al.lll996l ). 



'^Absorption cross section and electron affinity of circumovalene (EA=1.9 eV: lMalloci et al.ll2007l ) 
and maximum yield {Y = 1). 



is the FUV field measured in units of the iDraind ( 19781 ) field. Rate for Nc = 100. FUV and 
optical field using a T = 30, 000 K blackbody used for x > 1- 



qvan Dishoeck et al.l (|2006l ) except for the attenuation factors for H2O and OH photodissociation. 
Here we adopt Roberge et al (1991), see footnote h, because otherwise the models overproduce these 
molecules and O2 at the second peak compared with observations (see H09). 
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http://www.strw. leidenuniv.nl/ ewine/photo/index.php?file=pd.plip, Roberge et al (1991) 
'Rate for as s umed cross-section of 1 x 10^^^ cm^ and threshold of 957 A. Depth dependence from 



van Dishoeckl (119881 ) for a molecule with a photoionization threshold wavelength of 950 A. 



■^Ccrp is the primary cosmic-ray ionization rate per hydrogen nucleus. V arious rates are i nvest i- 
gated in this paper. The total rate including secondary ionizations is from iDalgarno et al.l (|l999l ). 

■^UDFAOe 



Brian k Mitchelll (|l99m 



McCall et al.l ^2004 ) 



"P. Goldsmith 2 011 private communication using mean total r ate from 



Jensen et al. ( 2000 ) of 5.7 x 10 ^ cm^ s ^; branching ratios from Jensen et al. ( 200ol ) 



et al.l torn and 



°Dielectronic plus radiative recombination rates from lBadnell et al.l ([2003) and iBadnelll (120061 ). 



POSU_01_2009 rate tables; http://www.physics.ohio-state.edu/~eric/research.html 



^State-specific rates from lStancil et al.l (jl999l ) 
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Table 2. Gas Phase Abundances and Grain Properties 



Species 


Symbol 


Value 


a 


Ref 


Carbon 


x(C) 


1.6 X 10" 


-4 


b 


Oxygen 


x{0) 


3.2 X 10" 


-4 


c 


Silicon 


x(Si) 


1.7 X 10- 


-6 


d 


Iron 


x(Fe) 


1.7 X 10- 


-7 


d 


Sulfur 


x(S) 


2.8 X 10" 


-5 


d 


Magnesium 


x(Mg) 


1.1 X 10" 


-6 


d 


PAHs 


x(PAH) 


2.0 X 10" 


-7 


e 


Grain Area 




2 X IQ-^^ 


f 



^Gas-phase abundances per hydrogen nucleus. 



Sofia et al.1 (|200J) 



Wever et al.l (|l998l ) 



nSavage Sz SembachI (jl996l ) cool diffuse cloud towards C 
Oph 



''Abundance from IWolfire et al.l (120031 ) modified for 
Nc = 100 planar PAHs. This abundance gives a total 
number of C in PAHs of 2 x 10^^ per hydrogen or ~ 6% 
of C in PAHs. 



^Units cm^ per hydrogen; iHollenbach et al.l ( 20091 ) 
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Appendix B: Simple Analytic Analysis of the Results 

The results in §3.1 and §3.2 can be understood by a simple analytic chemical model that 
incorporates the main physics. Such a model, though approximate, has the advantage of allowing 
one to determine and understand the sensitivity to various model parameters, and serves to validate 
the numerical model. 

Analytic solution to x(OH~^) and e in the first peak (Ay < Ij- The top panel of Figure 1 
describes the main chemical pathways to OH"*" for Ay < 1, and Table 1 lists the rate coefficients 
used for each of these reactions. Figure 1 omits the formation of OH^ by the photoionization of 
OH, which is only important at very low Ay ^ 0.01 and never dominates the production of the 
column of 0H+ for clouds of higher Ay- This minor route produces OH by the reactions O + 
H2 (vibrationally excited), O + H, and the formation of water ice on grains followed by desorption 
leading to OH. The photoionization of OH leads to a constant (low) OH"*" abundance at the cloud 
surface. We also neglect that route in the analytic solution, as it docs not affect the peak abundance 
of 0H+, nor the values of e near the peak. In addition, we approximate e for this peak by only 
including the rate of formation of 0H+ by the reaction H2 + 0+ — > 0H+ + H. 

In the following we use rate coefficients from Table 1 with two exceptions. The Stancil et al 
(1999) rate coefficients for the reaction of H+ with O and its reverse reaction are complicated, and we 
simplify them here, focussing on results for low density clouds. We set 71 = 4x 10~^°exp(— 230 K/T) 
cm^ as the rate coefficient for the reaction H"*" + O — > O"*" + H, and 72 = 4 x 10~^° cm^ s"^ 
as the rate coefficient for the reaction O"*" + H — ?> H"*" + O. Let 73 be the rate coefficient for the 
reaction H2 + O^ —?■ 0H+ + H; the rate coefficient for the recombination of H+ with electrons: 
ap- the rate coefficient for the neutralization of H"*" by PAH~; and ap the rate coefficient for the 
charge exchange of H+ with PAH. Let a;(H2) be the abundance of H2 with respect to H nuclei, Xg 
the abundance of electrons, x(PAH) the abundance of neutral PAHs, x(PAH~) the abundance of 
PAH~, and xq the gas phase elemental abundance of O. The solution to e, the efficiency of making 
0H+ from cosmic rays, is then: 

1 

where the parameter y is given: 



(1) 



^ [72 + 73a;(H2)] [aeXe + ap-g;(PAH-) + apx(PAR)] 

^ 11J3X0X(H.2) 

We utilized the fact that e approaches unity when y < 1 to set the condition in §3.1.1 (Eq. 1) on 
the abundances of O, H2, electrons, PAHs and PAH~s which lead to e ~ 1. One can immediately 
see that for low values of x(H2), where y > 1, e is proportional to a;(H2). Here, the dominant route 
of H+ destruction is by recombination with electrons, PAHs, and PAH~s. The small fraction of 
the H+ created by cosmic rays that charge exchanges with O to form 0+ that then reacts with 
H2 to form 0H+ is proportional to the H2 abundance. Substituting likely values of Xg ~ 3 x 10~^ 
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(wc focus on cases with high Ccrp/n-z and here as well as C"*" contribute electrons-see Figure 4), 
x(PAH)~ 1.85 X 10-^ x(PAH-)~ 1.5 x IQ-^, and xq ~ 3.2 x 10"^, we obtain 



e ~ 



227e-2-3/^2a;(H2) 



[1 + 4.25x(H,] [(3,^) T--- + 4.4 (f^) T,"- + 2.7 {^^) 



(3) 



where T2 = T/lOO K. For low values of x(H2) << 0.2, and standard values of Xg, 2;(PAH) and 
a;(PAH~), this becomes 



e ~ 28e-2-2/^2x(H2). (4) 

For higher values of x(H2) >> 0.2, e saturates at values near unity, as nearly every cosmic ray 
ionization leads to a production of 0H+. With standard values of Xg, a;(PAH) and a;(PAH~) 
Equation (3) then becomes 



^~ l + 0.15e2-3/T2- 

These formulae are very good fits to the results seen in Figures 4 and 6. Note, however, that they 
are only valid for x(H2) << 0.5, since we have implicitly assumed x(H)>> x(H2). 

Equation (3) shows that PAH~ in particular can lower e by neutralizing H+. We emphasize 
that equations (3-5) are valid only if the PAH and PAH~ rate coefficients we have adopted in Table 
1 are valid. Equations (4-5) assume the standard PAH and PAH~ abundances given above. 

Rather than a simple saturated value of e at high .'r(H2), however, we see in the figures that the 
first peak in 0H+ shows a value of e that peaks and then falls with increasing depth, once most of 
the gas is H2. This is because once the gas is more than half molecular (x(H2) > 0.25), the cosmic 
ray ionization is mostly of H2, which we have not included in the analytic treatment. The route to 
0H+ via the ionization of H2 is not as efficient at the surface, because of the high abundance of 
electrons which react quickly with H^. Thus, as x(H) declines deeper into the cloud, the value of 
e falls with increasing Ay- 

The analytic equation for a;(OH+) requires the addition of reactions of 0H+ with electrons 
and H2 that destroy 0H+. Let 74 be the rate coefficient for the destruction of 0H+ by H2, and 75 
be the coefficient for destruction by electrons. Then, we find 



x(OH+) = ^^-^^/""^ . (6) 

Here, Qrt ^ l-5Ccrp is the total rate of cosmic ray ionization, including secondaries. Ignoring the 
photodissociation term, we can again scale to likely values of parameters to obtain 
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x(OH+) 



3x10- 



j ) T^^-^ + 51x(H2) 



10-16 



crt \ I 100 cm" J 



s 1 / \ n 



Again, we can break this equation into two regimes. For a;(H2)<< 0.02(xe/3 x 10~^)T2^'^ , 0H+ is 
mainly destroyed by dissociative recombination with electrons and the second term in the denomi- 
nator can be ignored. Then, assuming standard parameters for the abundances of Xe, rE(PAH), and 
x(PAH~), we obtain 

.(0H+) 2.1 X [j^^) (^^^) • (8) 

However, if 0.02(xe/3 x 10-'^)r2^°-^ << x{B.2)< 0.2, then H2 mainly destroys 0H+ and, again using 
standard values for the parameters, we obtain the peak and plateau value of the OH^ abundance: 

x(OH-) . 4^1 X 10-.-/-. (^^) (15^) . (9) 

For higher values of x(H2) > 0.2, the abundance of 0H+ declines as e saturates, and as x(H) 
declines leading to the less efficient production of OH"'" initiated by cosmic ray ionization of H2. 
The column of OH"^ in the first peak can then be estimated 

iV(OH+) ^ 8 X lO^^e-^-^/^^ (^^^^) i^^v/^-^) (10) 

where AAy is the width of the peak region. 

This shows the important result that the peak abundance (and approximately the column of 
OH"'" in the first peak) is proportional to Ccrt/n. In order to find regions of high OH"'" column that 
can be observed, and where cosmic rays are the ultimate cause of the ions, we therefore must look 
at low density regions with high cosmic ray ionization rates. 

Analytic solution to the second peak of H^O^ in the regime Ay > 1 and with PAHs present. 
Deeper in the cloud {Ay > 1) the dust and gas shield the FUV photons and carbon is no longer 
in the form of but converts to CO. Here, the electron abundance drops rapidly with depth, 
and the gas is almost entirely molecular so that the route to HaO"*" is started with the cosmic ray 
ionization of H2. We then define e(H30"'") so that the rate of formation of H3O"'" per unit volume 
is Ccrtf^e(H30"^)El Again e(H30"'") is equivalent to the fraction of cosmic ray ionizations that lead 

to H3O"*". 

Figure 1 shows the route to H30"^ from the ionization of H2. The key competition that 
determines e(H30"'") revolves around Hj|". With all the hydrogen molecular, the routes via H2 



Note that if the total ionization rate per H is (,crt, then the total rate per H2 is roughly 2(^crt- The cosmic ray 
ionization rate per unit volume is then 2(^crt?i(H2) = C^crtn in fully molecular zones. 
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dominate the rates via electrons. Thus, the only interruption in the chain that leads to HsO"*" is 
the possibility that HJ^ will react with either electrons or CO rather than with O to form either 
0H+ or H20^. We then define 76 as the rate coefficient for reacting with O to form either 
OH"'" or H2O"'"; 77 is the rate coefficient for to dissociatively recombine with electrons; 78 is the 
rate coefficient for to react with CO; and 79 is the rate coefficient for H3O+ to dissociatively 
recombine with electrons (all channels). Therefore, 

e(H30+) c (11) 

^ ' ^ 76x(0)+77Xe+78s(CO) ^ ' 

The solution for the abundance of HaO"*" follows 

,(H30+) = (12) 

793^6 

Wc find in our model runs that electrons are formed by the cosmic ray ionization of H2 and 
destroyed by attachment to neutral PAHs. Most of the PAHs deep in the cloud are neutral, and 
therefore using a neutral PAH abundance of 2 x 10^^ we find 



Xe ^ 7.7 X 10-^(Ccrt/10"^*^ s~^)(100 cm-'-^/n) (13) 

If either O or CO dominate the destruction of H^ over electrons (or < 8 x 10~^a;(CO) + 6 x 
W-^x{0)), which we find is usually the case for most of our runs at the second peak of the HaO"*" 
abundance, then we find, substituting Xg into the equation for ^(HsO"*'), 

. 4 X 10-^7-/3(3 K)"-^^.r(0) 
^l^3U )- ^(o) + i.4x(C0) ^ ' 

At the peaJi, we generally find that most of the carbon is in CO, so x(CO) ~ 1.6 X 10-^ and that 
most of the remaining O is atomic, or x{0) ~ 1.4 x 10""^. Therefore, we predict a peak H3O+ 
abundance of 

a;p(H30+) ~ 1.6 x 10-*^(T/30 K)°-^^ (15) 

This explains why the second peak abundance of HsO"*" is independent of (^crt and n. The production 
rate per unit volume of HaO"*" depends on Ccrt^ but the destruction rate per unit volume depends 
on ri^Xe- Since electrons are formed by cosmic ray ionization and destroyed by PAHs, Xe oc (crt/n. 
Therefore, both production and destruction of HsO"*" is proportional to (crtn, and the resulting 
abundance of H3O+ is independent of both parameters. As one goes to higher Ay from the peak, the 
atomic O freezes out on grains as water ice, and thus the abundance of HaO"*" declines. The thickness 
of the region of peak H3O"*' abundance is therefore of order AAy ~ 1, or a hydrogen column of 
2 X 10^^ cm"^. Multiplying this by a;p(H30"^), we estimate columns of H3O+ of A^(H30"^) ~ 4 x 10^^ 
cm""^, as we found in our numerical runs. There has been some indication (Rimmer et al 2011) 
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that the cosmic ray ionization rate may decline with depth into a cloud. Prom the above result, 
this should have little consequence on the H3O+ column. 

We should emphasize that this prediction is dependent on the presence of PAHs deep {Ay ^ 5) 
in the cloud. If PAHs coagulate on grain surfaces, then PAHs do not remove free electrons from 
the gas phase. In this case, the electron abundance at the peak of the HsO"*" can be higher than 
in the PAH case if gas phase metal ions are present. This suppresses the HsO"*" abundance at the 
peak. Therefore, the observation of high columns of H3O"'" deep in molecular clouds may indicate 
either the presence of PAHs lowering the electron abundances there, or, in the absence of PAHs, 
may indicate a high depletion of gas phase metals at the second peak. 
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Appendix C: Sensitivity of results to other parameters 

Recently, Klippenstein et al (2010) have provided theoretical rate coefficients for the reactions 
of with O and with CO that are 10 to 20% higher than the rate coefficients we adopted. We 
have run our standard cases with these new coefficients and find that the columns of OH"*", H20''~, 
and HsO^ increase by roughly 10% with the new coefficients. In addition, Cartledge et al (2004) 
suggest a gas phase oxygen abundance of 2.84x10"^, whereas Jensen et al (2005) suggest a value 
of 4.2 to 4.7x10"^. We have adopted 3.2x10"^ but the variation in the literature suggests that we 
test the dependence of the ion columns on the gas phase abundance of elemental O. Assuming an 
abundance of 4.5x10"^, an increase of 40% over our standard rate, increases the columns of OH"*", 
H2O"'", and H3O+ by 10-20%. Therefore, the slight possible variations in these rate coefficients 
and/or the gas phase elemental O abundance have a negligible effect on our conclusions. 

We have also tested our results for their sensitivity to the H2 formation rate coefficient for 
formation on grain surfaces. Wc find that lowering this rate by a factor of 3 hardly affects the 
results for high density gas, and that the effect is most pronounced at low density. At n = 100 
cm~^, lowering the H2 formation rate by 3 lowers the columns of OH''" and 1120+ by factors of less 
than 2 in both peaks. The H3O+ column in both peaks decreases by a bit less than a factor of 
3. The main effect of lowering the H2 rate coefficient is to drive the H/H2 transition deeper into 
the cloud, which moves the first peaks of the ions to greater depths, but their columns and peak 
abundances do not change appreciably. 

In the text we discussed the effect of lowering the gas phase abundance of low ionization 
potential metal atoms such as S, Fe, Mg, and Si if PAHs are not present. Here, we examine the 
effect of raising their abundances. Our standard gas phase abundances for these species are given 
in Table 2 of Appendix A. In our test cases, we raise the gas phase abundances of these species 
to 10~^, a factor of 10-100 times their standard values. For our standard case of n = 10^ cm~^ 
and Ccrp= 2 x 10~^^ s~^, but with no PAHs, this raises the electron abundance at high Ay in the 
gas from about 10~^ to about 10~^. This rise depresses the (second) peak abundance of HsO"*" 
from about 10~^ to about 10"^'' if the photodesorption yields of Fe, Mg, and Si arc 10"^ a likely 
upper limit (H09). Smaller yields lower the electron abundance, thereby changing the HaO"*" less. 
However, if PAHs are present, metal ions recombine with PAH~ and PAH, which is a much faster 
process than with electrons. Hence, the electron abundance does not depend much on the metal 
abundances. Consequently, the (second) peak H3O+ abundance does not change significantly. The 
first peaks of the ions are not affected, because the electrons are supplied by C"*", whose gas phase 
abundance is more than 10 times the (high) abundance of the metal ions. 

Section 4 discusses observations that indicate A^(H30''") > 10^^ cm~^ in some clouds, possibly 
even as high as 10^^ cm^^. The only way our models can accommodate very large columns of HaO^ 
in a single cloud is if there is incomplete freeze-out of water ice on grains. Note that our models are 
steady state, so that time dependent effects might leave more elemental oxygen in the gas phase 
than steady state results would indicate. Another possibility is that the grains are sufficiently warm. 
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<: 100 K, to thermally evaporate water ice off the grains, or that cosmic ray rates are sufficiently 
high to desorb the ice mantles. We have run our code fixing the grain temperature to be > 110 K, 
so that water ice cannot form and elemental oxygen is plentiful in the gas phase at all Ay- In the 
case of PAHs, and for our standard case of n = 10"^ cm"^ and Ccrp= 2 x 10~^^ s~^, we find that 
x(H30+) plateaus at a value of about 2 x 10~^ at about Ay 6, and stays constant at that value 
for all higher Ay- (Note that given the discussion in §3.1, this result is independent of Qrp and n). 
Therefore, a cloud with a hydrogen nucleus column of 10^'^ cm~^ would have iV(H30"'")~ 2 x 10^^ 
cm~^, for example. Without PAHs, and assuming the low ionization potential metal atoms strongly 
deplete on grains, we would predict similar columns. 
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Fig. 1. — Standard ion-neutral chemistry leading to the formation of OH"*", H20~'', and H3O''' in 
clouds via ionization of atomic H (top panel), and via ionization of molecular H2 (bottom panel). 
In the top panel, the destruction of labelled "P/P~/e~" means the neutralization of H+ by 
reacting with PAHs, PAH~s, or electrons. "CR" means cosmic ray ionization. 
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n=102 cm-^ x=l. CrD=2xlO->^ s-> 




Fig. 2. — The variation of gas phase abundances of species as a function of depth Ay into a cloud 
for the standard case n = 100 cm~^, x = 1; and Ccrp= 2xl0~^^ s~^. To convert Ay to hydrogen 
nucleus column A^, use cm . This case probes diffuse cloud- like condition to high 

Ay, or could apply to low density surfaces of GMCs experiencing the local interstellar radiation 
field. Note these are columns when the cloud is viewed face-on. Color versions of the figures 
available in the on-line manuscript 
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n=102 cm-^ x=l. <'crD=2xlO->« s-> 




Fig. 3. — The variation of gas phase abundances of species as a function of depth Ay for the 
standard case n = 100 cm~^, x = Ij ^^'^ Ccrp= 2 x 10~^^ s~^. This is the same case as Figure 2 but 
with 10 times the cosmic ray ionization rate. This case probes diffuse cloud-hke condition to high 
Ay, or could apply to low density surfaces of GMCs experiencing the local interstellar radiation 
field. 
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Fig. 4. — The variation of the parameter e, the ratio of the rate of OH"*" formation to the cosmic 
ray ionization rate of H and H2, the gas temperature T (labehed on right), and the gas phase 
abundances of electrons, H2, Hj^, and OH"*" as a function of depth Ay for the case n = 100 cm~^, 
X = 1, and Ccrp= 2 x 10~^^ s^^. This case probes diffuse cloud-Uke condition to high Ay, or could 
apply to low density surfaces of GMCs experiencing the local interstellar radiation field. This is 
the same case as the previous figure. 
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n=10'' cm-3, x=10^ ^=2x10-'^ s-' 




Fig. 5. — The variation of gas phase abundances of species as a function of depth Ay into the 
cloud for the standard case n = 10^ cm~^, x = 100, and Ccrp= 2 x 10~^^ s~^. This case may be 
appropriate to GMCs with elevated FUV fluxes incident due to nearby O and B stars. 




Fig. 6. — The variation of the parameter e, the ratio of the rate of OH"*" formation to the cosmic 
ray ionization rate of H and H2, the gas temperature T (labehed on right), and the gas phase 
abundances of electrons, H2, , and OH"*" as a function of depth Ay for the case n = cm~^, 
X = 100, and Ccrp= 2 x 10^^^ s~^. This is the same case as the previous figure. 
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n=10« cm-3. x=10^ ^=2x10-^^ s-' 




Fig. 7. — The variation of gas phase abundances of species as a function of depth Ay for the case 
n = 10^ cm~^, X = 10^) ^iid Ccrp= 2xl0~^^ s~^. This case demonstrates the structure of PDRs 
with both high density and high FUV fluxes, where elevated (> 300 K) temperatures in the region 
with significant H2 leads to the enhanced production of H"'" by chemical routes not initiated by 
cosmic ray ionization (see text). One mark of this is the enhanced OH abundance, produced by 
the neutral- neutral reaction of H2 with O that is seen at Ay ~ 1, where the gas temperature is 
T ~ 1000 K. The enhanced OH reacts with to produce CO"'', which then reacts with H to form 
H+. The rest of the chemistry leading to 0H+, H2O''', and HsO"*" is seen in Figure 1. 
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Iog,o[N(OH-)/cm-2]. Cp=2x10->« s" 




logio[Vcm-3] 

Fig. 8. — Contours of the column of OH^, labelled in log units, as a function of n and x for a fixed 
primary cosmic ray ionization rate of Ccrp=2 x 10~^^ per H atom. The upper left portion of 
the figure is blacked out because radiation pressure on dust drives dust quickly through the PDR 
in this region, invalidating the physics assumed in the model. This combination of x ^'^d n are 
rarely observed in any case. Except in the upper right hand corner of this figure (high n and high 
x), we see that for fixed cosmic ray ionization rate, the column is roughly proportional to n"^, 
and independent of %. The upper right hand corner shows a secondary peak in the 0H+ column, 
caused by the alternate chemical routes described in text and shown in more detail in Figure 7. 
Low densities n < 100 cm~^ are required to obtain columns greater than about 10^^ cm~^ created 
by cosmic rays. 
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logio[Vcm-3] 

Fig. 9. — Contours of the column of H20^, labelled in log units, as a function of n and x foi' ^ 
fixed primary cosmic ray ionization rate of Ccrp=2 x 10~^^ s~^ per H atom. The same discussion 
as in the previous figure applies here. To obtain columns greater than about 10^^ cm~^ created by 
cosmic rays requires low densities n < 100 cm^^. 




Fig. 10. — The columns of OH"*", A^(OH^), are plotted as a function of Ayt for four values of the 
primary cosmic ray ionization rate per H atom divided by n2 = n/100 cm~^ and for x/"-2 = 3.16. 
The dotted line plots the abundance of H2 at cloud center, and these values appear on the right of 
the figure. Recall Ayt is the total Ay through the diffuse or translucent cloud. The H2 abundance 
at cloud center, x(H2), is also plotted (dotted lines) and its values noted on the right. 
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Fig. 11. — The columns of OH"*', A^(OH''"), are plotted as a function of Ay^ for four values of the 
primary cosmic ray ionization rate per H atom divided by n2 = n/100 cm~^ and for x/''^2 = 10. 
The dotted line plots the abundance of H2 at cloud center, and these values appear on the right of 
the figure. This figure is the same as Figure 10, only with raised by 3.16. 
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Fig. 12. — The ratio A^(0H^)/A^(H20''') is plotted as a function of Ayt for four values of the 
primary cosmic ray ionization rate per H atom divided by = n/100 cm~^ and for x/"-2 = 3.16. 
The dotted line plots the abundance of H2 at cloud center, and these values appear on the right of 
the figure. 
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Fig. 13. — The ratio A^(0H^)/A^(H20''') is plotted as a function of Ayt for four values of the 
primary cosmic ray ionization rate per H atom divided by n2 = n/100 cm~^ and for x/''^2 = 10. 
The dotted line plots the fraction of H2 in the cloud, and these values appear on the right of the 
figure. This figure is the same as Figure 12, only with x/''^2 raised by 3.16. 
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x/n2=3.16 




N(0H+)/N(H20+) 

Fig. 14.— Logio [N{OR+)/[J (n(H)/T2) dz] is plotted on the vertical axis and iV(0H+)/iV(H20+) 
on the horizontal axis for the case x/'^2 = 3.16. Plotted as solid lines are constant values of 
logio[^yt], labelled on the bottom of these lines. Plotted as dashed lines are contours of constant 
logio[Ccr-p/'^2], labelled on the left, and in units of s"^. The two data points at right are two velocity 
components (diffuse clouds) toward W49N. The lower limit data point is toward W31C (see text). 
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Fig. 15. — This figure is identical to Figure 14 except that x/"2 = 10. 
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X/n2=3.16. no PAH/PAH" + H* 




N(0H+)/N(H20+) 

Fig. 16. — This figure is identical to Figure 14 = 3.16) except that the PAH and PAH— rates 

with H^ have been reduced by ^ 4 so that H"*" is mainly destroyed by electrons or by forming O"*" 
which then reacts with H2 to form 0H+. Note that the inferred Ccrp/'^2 values decrease compared 
to the case with PAHs. 



-61 - 




Fig. 17. — This figure is identical to Figure 15 = 10) except that the PAH and PAH— rates 

with H^ have been reduced by ^ 4 so that H"*" is mainly destroyed by electrons or by forming O"*" 
which then reacts with H2 to form 0H+. Note that the inferred Ccrp/'^2 values decrease compared 
to the case with PAHs. 



